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PREFACE 

This  report  discusses  a  slxople  model  of  long-range  sound  trsas- 
alsslon.  It  Is  part  of  RAND' s  continuing  study  of  the  atmospheric 
waves  generated  by  nuclear  explosions.  A  later  report  Is  planned 
which  will  apply  the  theory  developed  In  this  report  to  the  e:q>lana- 
tlon  of  some  of  the  characteristic  features  of  the  waveforms  detected 
at  large  distances. 

This  report  should  be  of  inteirest  to  agencies  and  contractors 
concerned  with  the  detection  of  nuclear  esqploslons.  It  should  also 
be  of  Interest  to  other  reseeurch  wozicers  who  are  studying  the 
phenomena  associated  with  wave  propagation  In  Inhomogeneous  media. 
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SUWAHY 

A  theoretical  study  has  been  made  of  the  transient  sovind  propa¬ 
gation  from  a  point  source  In  a  triple-layered  medium  consisting  of  a 
homogeneous  fluid  layer  saxidwlched  between  two  similar  half  spaces  of 
the  same  density  but  of  higher  sound  speed,  nihe  source  spectrum 
Is  such  that  a  large  number  of  normal  modes  sure  excited.  Both  source 
and  receiver  are  assumed  to  be  In  the  same  half  space.  The  generc^. 
solution  consists  of  a  direct  wave,  a  wave  reflected  from  the  near 
edge  of  the  center  layer,  a  lateral  wave,  and  the  normal  and  complex 
mode  waves.  At  large  distances  r  the  normal  mode  waves  dominate. 

It  is  fovind  that  one  may  relate  frequency,  pheise  velocity,  group 
velocity,  and  excitation  amplitudes  parametrically,  and,  on  the 
basis  of  this,  the  variation  of  the  characteristics  of  tie  normal 
modes  with  mode  number  n  is  studied,  similar  methods  can  be  used 
to  study  the  conqtlex  modes.  Each  normal  mode  is  found  to  carry  a 
limited  band  of  frequencies — the  beuad  width  for  higher  order  modes 
being  8q>proxlmately  proportional  to  n~^.  Unless  r  is  e:rtremely 
large,  the  amplitudes  of  the  Airy  phases  for  the  higher  order  modes 
will  be  negligible.  On  the  basis  of  the  study,  one  may  give  a 
qu8Q.ltatlve  description  of  the  waveforms.  One  effect  of  placing 
source  and  receiver  outside  the  center  layer  Is  to  accentuate  lower 
order  modes  and  the  lower  frequency  portion  of  the  source  spectrum. 
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I.  HITRODUCTIOII 

In  the  hierarchy  of  guided  wave  problema  In  stratified  media, 
that  of  giving  a  theoretical  interpretation  to  the  transient  vave- 
foxms  observed  at  large  dlstcmces  from  a  high  frequency  source  Is 
generally  considered  to  be  among  the  more  difficult.  Although 
both  of  the  stemdard  mathematical  approaches  generally  utilized 
In  guided  wave  problons  (namely  the  so-called  nomal  mode^^~^^  and 
ray  theory  approaches may  in  principle  be  applied,  neither 
seems  to  lead  to  an  easily  Interpretable  solution.  Ihe  reasons 
for  this  difficulty  are  readily  understandable.  Since  the  distance 
frcan  the  source  Is  large,  the  number  of  Individual  rays  which  must 
be  summed  will  also  be  large.  Also,  since  the  frequency  band  of 
the  source  lies  above  the  cutoff  frequencies  of  a  large  niniber  of 
normal  modes,  one  must  generally  siqperimpose  many  modes. 

Intuitively,  one  would  feel  that  the  normal  mode  approach  is 
the  better  of  the  two,  as  It  Is  an  asymptotic  theory  which  enables 
one  to  separate  out  at  the  very  beginning  that  portion  of  the  wave 
which  will  predominate  at  large  distances.  However,  for  even 
relatively  simple  models  of  stratified  media,  the  calculation  of 
the  pertinent  characteristics  (l.e.,  group  velocities  and  excita¬ 
tion  amplitudes)  of  the  normal  modes  must  be  done  nunerically. 

(ihe  only  exception  of  ^Ich  we  are  aweo^  Is  the  model  In  which  the 
sound  speed  varies  as  in  HPBteln's  symmetrical  layer. As 
long  as  one  Is  Justified  In  considering  one  or  a  small  number  of 
modes,  this  does  not  present  insurmountable  difficulties.  Indeed, 
some  very  convenient^  methods  have  been  devised  for  accosqpllshlng 
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this  for  models  of  various  sophistication.  Nevertheless >  within 
this  context >  the  study  of  the  systematic  variation  of  mode  character¬ 
istics  with  mode  number  requires  the  numerical  calculation  of  the 
characteristics  for  a  lai^e  number  of  modes. 

To  minimize  such  difficulties,  It  Is  convenient  to  begin  with 
the  study  of  a  very  sloqile  model  and  to  consider  this  model  In  some 
detail.  The  purpose  of  this  paper  Is  to  report  on  the  results  of 
such  a  study. 

The  model  ve  chose  was  that  of  a  homogeneous  fluid  layer 
sandwiched  between  two  similar  homogeneous  half  spaces  of  higher 
sound  velocity.  Mathematically,  this  model  Is  very  slmlleur  to  that 
originally  proposed  by  Pekerls^^^  for  the  study  of  sound  propagation 
In  shallow  water.  For  example,  the  dispersion  characteristics  of 
the  second,  foujrth,  etc.,  modes  for  our  model  are  Identical  to  those 
of  the  first,  second,  etc.,  modes  for  Pekeris*  model.  Our  principal 
reason  for  considering  this  particular  model  rather  than  Pekerls' 

Is  that  our  Interest  Is  primarily  in  sound  propagation  In  the 
atmosphere  rather  than  In  the  ocean,  and  we  felt  that  our  model 
would  hold  a  greater  relevancy  to  this  subject. 

Whether  or  not  this  model  Is  relevant  to  the  understanding  of 
sovind  propagation  In  the  atmosphere  Is  questionable  at  the  present. 

Our  belief  In  Its  relevancy  Is  based  on  the  fact  that  the  first 
temperature  mlnlnnim  In  the  atmosphere^ lies  some  distance  20  km) 
above  the  ground,  and  that  the  temperature  maximum  at  the  ^  km 
inversion  Is  less  than  that  at  the  groiind.  It  therefore  seems 
reasonable  to  believe  that  higher  frequency  waves  trapped  In  the 
lower  atmosidierlc  sound  channel  will  be  little  affected  by  tie 


-3 


presence  of  the  ground.  Naturally,  this  will  not  be  true  for  the 

very  lew  frequency  waves  (the  so-eslled  gravity  waves),  which  have 

C 14-17) 

recently  been  studied  by  a  nunber  of  workers,^  '  but  should 
become  of  Increasing  validity  when  one  considers  higher  frequencies. 

A  model  such  as  ours,  though  vastly  overs Isgpllf led,  might  be  help¬ 
ful  In  acquiring  a  qualitative  understanding  of  some  aspects  of 
the  physical  phenomena.  In  the  present  paper,  however,  we  shall 
not  e]q>lore  the  possibility  of  the  relevancy  of  our  model  to  this 
subject . 

To  simplify  our  model  as  much  as  possible,  we  have  assmned  that 
the  middle  layer  and  surrounding  half  spaces  etll  have  the  same 
density,  though  we  have  otherwise  (to  as  large  an  extent  as 
practicable)  avoided  any  numerical  specification  of  the  ratio  of 
the  two  sound  speeds  or  the  width  of  the  center  layer.  We  have  tried 
to  treat  the  problem  in  sane  generedlty  and  to  make  maxlsium  \ise 
of  analytic  techniqxies  rather  than  to  rely  solely  on  numerical 
calculations.  For  the  most  part,  we  have  limited  our  consideration 
to  the  case  where  both  source  and  receiver  are  outside  the  channel 
and  on  the  same  side  of  the  channel  (l.e.,  where  both  lie  In  the 
same  half  space).  The  mathematical  formulas  obtained  for  other 
source  and  receiver  posltlonlngs  are  appreciably  different,  and  It 
was  felt  that  a  coxtslderatlon  of  all  possible  cases  wotild  unduly 
lengthen  the  paper. 

This  specification  was  also  motivated  by  the  desire  that  o\ir 
model  have  some  relevancy  to  the  study  of  propagation  In  the 
atmosphere— where  both  source  and  receiver  are  commonly  near  the 


ground  and  hence  (In  a  certain  senee)  outside  the  channel.  This 
departs  from  the  usual  procedure  of  placing  both  Inside  the  channel 
as  Is  done  In  the  majority  of  paqpers  on  wave  propagation  In  the 
ocean.  The  configuration  ve  consider  has  received  little  attention 
although  Officer  does  briefly  discuss  this  tqplc  In  his  book.^^^^ 

In  our  analysis  of  the  problem,. ve  have  been  aided  by  the  fact 
that  the  characteristics  of  the  normal  modes  may  all  be  related 
parametrically  to  the  frequency  by  the  Introduction  of  a  variable 
U.  This  Is  possible  even  If  one  takes  the  density  of  the  center 
layer  to  be  different  from  that  of  the  surrounding  half  spaces, 
although  considerable  slnqpliflcatlon  Is  achieved  by  taking  all 
densities  equal.  Studying  the  parametric  equations  enables  us  to 
derive  several  properties  concerning  the  mode  characteristics  and 
furnishes  us  with  a  convenient  method  of  studying  the  variation  of 
these  characteristics  with  mode  nuinber. 

Placing  the  source  and  receiver  outside  the  channel  turns  out 
to  have  some  Interesting  consequences.  The  excitation  aoplltudes 
aie  of  magnitude  over  only  a  limited  range  of  fre¬ 

quencies  for  each  mode,  and  the  overlap  of  the  frequency  bands  for 
different  modes  Is  small.  For  each  normal  mode  the  snplltudes  of 
the  "ground  vave"  are  generally  much  larger  than  those  of  the 
"water  vave"  or  the  Airy  phase  unless  the  distance  from  the  source 
to  receiver  Is  exceptionally  large.  The  maximum  amplitude  of  each 
mode  decreases  with  Increasing  mode  number,  giving  the  Impression 
that  the  model  favors  lower  frequencies.  These  effects  shall  be 
discussed  In  some  detail. 
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II.  THE  POHMAL  90IAJTK« 

Tba  model  ve  shall  consider  Is  eonv«niently  described  by 

specifying  the  dependence  of  the  speed  of  sound  c  on  the  vertical 

coordinate  *.  For  z  >  h/2  or  z  <  -  H/2  (i-e.,  In  the  two 

half  spaces)  the  sound  speed  Is  c^.  In  the  center  layer,  where 

-  h/2  <  z  <  1^2,  the  sound  speed  Is  Cj^>  where  c^^  ^  ®o* 

density  p  Is  constant  throughout  the  media.  A  point  source  of 

sound  Is  located  on  the  z-axls  at  a  distance  z  from  the  plane  of 

o 

symmetry  (Fig.  1).  In  the  Inmedlate  vicinity  of  the  source  the 
velocity  potential  behaves  like  an  outgoing  spherical  wave, 

t  »  f(t-R/c)  /r,  (2.1) 

where  R  Is  the  distance  from  the  source  and  f(t)  Is  the  time  de¬ 
pendence  of  the  source.  In  general,  Is  the  solution  of  the  simple 
wave  equation  and  Is  continuous  with  continuous  gradients  everywhere 
except  at  the  source. 

The  techniques  of  finding  a  formal  solution  to  problems  of 

(1-4) 

this  type  have  been  discussed  at  some  length  in  the  literature.' 

The  general  procedure  Is  to  ei^ress  It  as  a  double  Integral  over 
frequency  u  and  wave  number  k.  In  treating  the  case  when  both 
soiirce  and  receiver  are  In  the  same  hsdf  space,  however.  It  Is 
convenient  first  to  extract  a  direct  wave  and  a  wave  reflected 
from  the  near  edge  of  the  guide,  such  that 
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♦  -  f(t-(Rp/c^,))/Rjj 
-  f(t-(Rj/c^))/Rj 


+  i,> 


(2.2) 


vbsre 

„  r  2  /  x2,l/2 

Rjj  =  [r  +  (z-z^)  1  ' 
and 

Rj  =  [r®  +  (Iz+z^l  - 

The  extraction  conforms  to  our  Intuitive  notion  that  these  waves 

should  he  present  and  aids  In  shortening  our  presented  Ion.  A 

tentative  Justification  for  assigning  a  negative  sign  to  the  second 

term  can  he  given  if  one  notes  that  the  Rayleigh  reflection  co- 
(19) 

efficient  for  angles  near  grazing  incidence  Is  nearly  -1.  Since 
we  are  Interested  In  the  case  of  large  r,  any  wave  reflected  from 
the  near  edge  of  the  guide  must  he  considered  u  being  near  grazing 
Incidence. 

After  the  extraction  of  the  two  waves,  the  remaining  term  <ir' 
may  he  written 

*•  =  2Re  J  g(u)e"^‘**du  T  (kr)Ukdk,  (2.5) 

*  "CO 

O 

where  g(u)  Is  the  Fourier  transform  of  f(t),  and  U  Is  independent 
of  the  radial  coordinate  r.  For  our  assumed  receiver-source  con¬ 
figuration,  U  Is  essentially  a  function  of  three  quantities: 
and  0,  where 
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(2.4a) 

4  =  -  k®, 

(2.4b) 

D  »  1  (!a+Zp|  -  H). 

(2.5) 

Tbe  phase  of  (or  is  taken  to  be  0  or  ^  depending  on  whether 
k  Is  less  or  greater  than  (w/c^)  (or  (u/cj^)  ).  The  quemtlty  D 
represents  the  sum  of  the  distances  of  the  source  and  receiver 
from  the  near  edge  of  the  guide  In  units  of  l/2  H*  Scpllcltly^ 


U 


1  (1-va^) 

ii-yTcr) 


i3pDH/2 


(2.6) 


where 


V  =  (Po  -  *  Pi) 

(2.7) 

Q  =  e 

(2.8) 

The  derivation  of  (2.6)  proceeds  along  standard  llnes^^~^^  and  Is 
omitted  In  the  Interests  of  brevity. 

The  Integration  over  k  In  Eq.  (2. 3)  Is  properly  a  contour 
Integration  In  the  complex  plane.  To  Insure  causality^  the  path 
must  pass  above  eOl  poles  lying  on  the  negative  real  axis  and  be¬ 
low  all  poles  lying  on  the  positive  real  axis.  Similarly,  the 
Integration  over  u  In  (20)  is  understood  to  pass  above  all  poles 
on  the  real  axis. 

Since  the  normal  mode  method  rests  on  a  contour  defoxmatlon 
In  the  k  plane,  It  Is  necessary  to  extend  the  definition  of  the 


-8- 


Integrand  to  Include  coqplex  values  of  k.  Xhls  Is  aeeoapllshed 

by  placing  branch  lines  and  requiring  that  the  Integrand  be  analytic 

everyvhere  except  at  these  branch  lines  and  at  Its  poles.  The 

placing  of  the  branch  lines  Is  standard.  One  branch  line,  made 

necessary  by  the  definition  of  the  Hankel  function,  extends 

vertically  dovnvards  from  the  origin.  Those  necessaiy  for  the  unique 

specification  of  the  phase  of  0^  extend  vertically  upwards  from  the 

branch  point  at  k  =  u/c^  and  vertically  downwards  from  k  =  -  u/c^. 

Branch  lines  for  0^are  not  essential,  as  U  Is  Invariant  If  0^  Is 

replaced  by  -  0^.  However,  since  U  Is  not  explicitly  given  as  a 
2 

function  of  0^,  It  Is  convenient  to  add  branch  lines  extending 
upwards  from  u/c^  and  doimwards  from  -  w/c^  to  avoid  any  inconsisten¬ 
cies.  It  Is  readily  seen  that  our  previous  specification  of  the 
phases  of  0^  and  0^  along  the  real  axis  Is  consistent  with  the 
placing  of  our  branch  lines. 

With  the  contoxir  deformation  appropriate  to  the  method  of 
normal  modes,  the  path  In  the  k  plane  Is  transformed  Into  a  sum  of 
contour  Integrals  around  the  poles  lying  In  the  upper  half  plane 
and  on  the  positive  real  cucls  plus  a  branch  line  Integral  around 
the  line  extending  upwards  from  u/c^.  This  contour  is  sketched 
In  Fig.  2. 

With  this  deformation,  the  quantity  j  may  be  decomposed  Into 
three  terms, 

*  *0,  Ml  -  (2-9) 

where  the  first  term  (normal  mode  solution)  represents  the  contribu¬ 
tion  from  the  poles  on  the  real  axis,  the  second  term  (complex  mode 


solution)  represents  the  contribution  from  the  poles  In  the  upper 
half  plane,  and  the  third  term  (lateral  wave)  represents  the  con¬ 
tribution  from  the  branch  line  IntegreiL. 

The  contovir  Integrals  around  the  poles  may  be  formally  evalua¬ 
ted  by  the  method  of  residues.  A  simplified  expression  for  the 
residues  may  be  derived  If  one  notes  that  the  poles  coincide  with 
the  roots  of  the  equation 

1  -  =  0.  (2.10) 
Then,  one  may  show  that,  near  a  pole,  U  Is  of  the  form 

^  “  kH(k-kp)  ' 

where  k_  denotes  the  location  of  the  pole  and 
P 


(2.12) 


Is  to  be  evaluated  at  k  =  k^.  !nie  result  of  the  Integration  around 
a  pole  may  readily  be  obtained  by  taking  2x1  times  the  limit  of 
the  product  of  (h-k^)  and  the  Integrand  as  k  approaches  k^. 

Since  the  positions  of  the  poles  will  generally  change  con¬ 
tinuously  with  u,  each  pole  may  be  assigned  an  Index  n  (for  those 
on  the  real  axis)  or  m  (for  those  above  the  real  axis),  where  n, 
m  =  1,2,?,  etc.  The  manner  of  assigning  these  Indices,  or  mode 
nunibers,  may  be  made  unambiguous  and  will  be  described  In  the  next 
section.  As  is  well-known,  poles  corresponding  to  a  given  mode 
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nuniber  vlll  not  necessarily  be  present  for  all  values  of  u.  in 

general  they  vlll  be  present  over  some  range  of  frequencies^  and 

one  may  apeak  of  a  lov  frequency  or  high  frequency  cutoff  (u  or 

c 

up  for  each  mode.  The  poles  on  the  real  axis  (l<e.,  the  normal 
modes)  vlll  have  no  high  frequency  cutoff,  but  those  off  the  real 
axis  (l*e.,  the  conqplex  modes)  vlll  In  general  have  both  lov  fre¬ 
quency  and  high  frequency  cutoffs.  (Die  latter  shall  be  shovn  to 
be  true  In  the  next  section  and  In  Appendix  A.) 

Up  to  nov,  the  method  ve  have  Indicated  for  evaluating  the 
normal  mode  or  coiiq)lex  mode  solutions  Is  first  to  perform  a  sum 
over  poles  and  then  Integrate  over  frequency,  with  the  Introduc¬ 
tion  of  the  concept  of  mode  nuiift>er  It  Is  possible  to  Interchange 
this  sequence.  Thus  the  formal  expression  for  may  be  given 
In  the  fom 


4jc  * 

♦hM  “  H"  ^  I  J  10H^^^(kr)g(u)e"^'^du,  (2.13) 

n^l  w 
c 

vhere  k  Is  h^(<<>)  and  the  cutoff  frequency  u^  depends  on  n.  A 

similar  expression  holds  for  the  chief  formal  distinction  being 

that  the  Integration  for  each  mode  has  a  high  frequency  cutoff  u' . 

c 

The  remaining  portion  of  the  solution,  the  lateral  vave.  Is 

expressible  In  the  form 

^  u/c^+l* 

♦l  “  J  S(u)e‘^‘*du  J  [U(p^)-U(-p^)]H^^^(kr)kdk, 

0  u/c 

'  o 


11 


where  we  have  chcmged  the  integration  over  k  from  the  original  path 
(which  proceeded  down  one  side  of  the  brcuoch  line  and  up  the  other) 
to  a  path  from  the  branch  point  upwards  along  the  right  side  of  the 
branch  line.  This  was  accomplished  by  making  use  of  the  fact  that 
the  only  formal  distinction  in  the  Integrand  on  the  two  sides 
of  the  branch  line  was  a  difference  of  it  in  the  phase  of  The 
value  of  to  be  utilized  in  (2.lU)  is  that  on  the  right  side 
of  the  branch  line. 

The  decomposition  of  the  velocity  potential  into  a  direct  wave, 
reflected  wave,  normal  mode  solution,  complex  mode  solution,  and 
lateral  wave  as  outlined  above  is  (subject  to  some  mathematical 
fine-points)  a  rigorous  expansion  of  the  solution  along  the  lines 
of  the  normal  mode  approach. 
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III.  lOCATIOH  OF  THB  POLES  IM  THE  COMPLEX  FLAME 

One  'beurrler  In  the  utilization  of  the  formal  solution  given 
in  the  previous  section  is  that  of  finding  the  expressions 
l.e.^  of  locating  the  positions  of  the  poles.  These  are  the  roots 
of  the  transcendental  equation  (2.lo)>  for  which,  in  general, 
explicit  solutions  in  terms  of  elemental^  functions  cannot  be  ob¬ 
tained.  Ve  shall  here  introduce  a  method  which  will  facilitate 
the  study  of  these  roots  and  which  is  convenient  for  numeric  clL 
computations.  We  shall  not  restrict  our  attention  to  the  real  roots 
alone  (although  these  are  of  primary  Interest),  as  our  method  makes 
a  study  of  the  complex  roots  feasible  without  an  undue  amo\ut  of 
effort . 

2  2 

Let  us  first  note  that  the  difference  is  Independent 

of  the  wave  number  k.  In  partlculcu:,  we  have 

^1  ■  ®o  =  =  («/c^)V,  (5.1) 

where 

5  -  (2/h)c^c^(c^  -  c®)"^/®  (5.2) 

is  a  frequency  and 

is  a  dimensionless  constant  which  characterize  the  model.  (Both  A 
and  u  are  real  since  we  are  assuming  ^  ^2.'^ 
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The  lack  of  dependence  of  (5.I)  on  k  euggcsta  the  follovlnf; 
svkbstltutlon: 

=  (2b/h)  cos  (1  (3.*«i) 

=  i(2B/H)  sin  n  ,  (3«^t) 

where  we  have  abbreviated  B  »  u/u.  In  terns  of  the  variable  m  Just 
Introduced,  the  equation  (2.10)  would  appear  ets 

exp  [41(3  cos  M  -  ^l)]  =  1,  (3*5) 

which  has  an  appealing  sls^liclty  when  contrasted  with  Its  original 
form  In  terms  of  k. 

There  cu:e  certain  mathematical  pitfalls  Involved  In  the  utili¬ 
zation  of  (3.3)  Instead  of  (2.10),  and  some  care  should  be  exer¬ 
cised.  We  shall  first  restrict  the  real  part  of  w  to  lie  between 
-ff  and  X.  With  this  restriction  and  the  definitions  (3*4) >  a 
unique  value  of  u  Is  associated  with  each  point  In  the  k-plane. 

We  shall  consider  this  association  In  some  detail. 

The  k-plane  may  be  Imagined  as  consisting  of  eight  regions — 
the  boundaries  of  these  regions  being  formed  by  the  real  and 
Imaginary  axes  and  the  branch  lines  extending  from  +  u/c^  and  +  u/c^. 
Although  It  is  laborious  to  make  a  detailed  specification  of 
and  p^  for  every  point,  it  is  a  relatively  simple  matter  to  specify 
the  quadrants  In  which  their  phases  lie  for  each  of  the  eight  regions. 
For  example,  in  the  region  to  the  right  of  the  branch  line  extend¬ 
ing  from  u/c^  and  above  the  real  axis,  the  phases  of  both  p^  and 
p^  will  lie  In  the  second  quadrant.  Providing  this  specification 


Is  a  routlns  problem  In  the  theory  of  conplex  variables.  Further¬ 
more,  If  ve  know  the  q]uadrants  In  which  these  phases  lie,  we  may 
determine  via  equations  (3*^)  the  sign  of  the  Imaginary  part  of  u 
and  the  quadrant  In  which  Its  reed  part  lies.  In  the  region  mentioned 
above,  for  exanple,  Is  negative  and  Is  In  the  second  quadrant. 

We  have  carried  out  this  analysis  for  each  of  the  eight  regions. 

The  results  are  given  In  Fig.  3. 

Let  us  now  note  that  equation  (3.3)  may  be  satisfied  only  If 
the  Imaginary  part  of  6  cos  u  -  |i  Is  zero.  When  one  writes  out  the 
Imaginary  part  e^llcltly,  he  discovers  this  can  be  true  only  If 
either  *  0  or  sin  Is  negative.  Referring  to  Fig.  3>  we  see 
that,  if  a  root  of  (3.5)  is  to  correspond  to  a  value  of  k  which 
satisfies  (2.10),  this  value  of  k  must  lie  either  on  the  real  axis 
between  «/ c^  and  u/ c,  or  between  -«/ c^  and  -«/ c, ,  or  else  it  must 
lie  In  one  of  the  two  regions  which  are  shaded  In  Fig.  3.  ISius 
when  studying  the  real  roots  of  (2.10),  we  shall  restrict  our 

attention  to  roots  of  (3'3)  for  which  u-  =  0  and  u  Is  between  0 

X  R 

and  x/2.  when  studying  the  conqplex  roots  we  restrict  our  attention 
to  roots  of  (3.5)  for  which  (x_<  0  and  Is  between  -  x/2  and  0. 

since  (2.10)  Is  Invariant  under  replacing  k  by  -k,  the  roots 
are  synmetrlcally  located  In  the  k-plane.  We  shall  accordingly 
limit  our  consideration  to  those  roots  In  the  first  quadrant  and 
on  the  positive  real  axis.  By  utilization  of  either  (3.4a)  or 

O  p 

(3.4b)  and  the  definitions  of  and  one  may  formally  ejqpress 
2 

k  In  terms  of  m: 


k®  -  («/Cp)®(l  +  A®  sln^u). 


(3.6) 
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It  |Ji  lies  In  either  of  the  two  regions  defined  In  the  preceding 
porsgrsphj  the  shove  equation  will  lead  to  a  unique  value  of  k  In 
the  first  quadrant  or  on  the  positive  reed  axis.  If  u  Is  reed  and 
between  0  and  n/2.,  the  value  of  k  computed  In  this  manner  will 
always  lie  between  u/c^  and  “/c^,  euad  any  such  root  of  (3*5)  will 
correspond  to  a  root  of  (2.10).  On  the  other  hand.  If  0  and 
-  it/2  <  <  0>  the  k  computed  from  (5.6)  will  not  necesseurlly  lie 

to  the  left  of  the  branch  line  which  extends  from  u/c^.  If  it  lies 
to  the  right  of  the  branch  line,  the  equations  (3*^)  will  be  in¬ 
compatible  with  the  definitions  of  and  This  possibility  can 
occur,  since  every  ^  will  not  necesseurily  correspond  to  a  point 
in  the  k-plane. 

The  correct  procedure  for  utilizing  (3-5)  is  readily  apparent. 

To  find  all  of  the  roots  of  (2.lo)  lying  on  the  positive  real  axis 

one  first  finds  all  the  roots  of  (3*5)  for  which  u  is  real  and 

between  0  and  a/2,  and  then  computes  k  from  (3.6).  To  find  all  of 

the  complex  roots  lying  in  the  first  quadrant  one  first  finds  all 

roots  of  (3*5)  for  which  u_<  0  and  -a/2  <  u  <  0  and  computes  a  k 

X  R 

lying  intite  first  quadreuat  via  (3>6)  for  each  such  He  then 
eliminates  all  complex  k's  which  do  not  lie  to  the  left  of  the 
branch  line.  It  la  easily  verified  that  the  roots  obtained  in  this 
manner  and  their  negatives  (inversions  through  the  origin)  constitute 
all  of  the  roots  of  (2.10). 

It  should  be  mentioned  at  this  point  that  the  parameter  has 
no  immediately  apparent  physical  significance.  It  is  merely  a 
parameter  which  aids  in  studying  the  solution  of  the  problem. 
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rae  Beal  Roota 

Bq>Mitlon  (3<5)  Is  satisifled  only  If  B  cos  -  4  Is  a  multiple 
of  x/2.  Thus  4  must  satisfy  an  equation  of  the  form 

B  cos  u  -  ^  =  (L/2)(n-l)jt,  (3*7) 

where  n  Is  an  Integer.  Ve  shall  consider  to  be  real  and  shall 
tentatively  Identify  n  as  the  mode  number  for  the  normal  modes. 

There  will  be  solutions  of  (3 >7)  between  0  and  ](/2  only  If 
n  >  1  and  B  >  l/2  (n-1).  The  low  frequency  cutoff  for  the  n-th 
normal  mode  Is  therefore 

“c  °  I  {n-l)*ij.  (3*8) 

If  n  Is  cm  even  Integer,  this  will  be  the  same  as  that  for  the 
n/2-th  mode  In  Fekerls*  model. 

If  u  lies  between  the  cutoff  frequencies  of  the  n-th  and 
(n  +  l)-th  modes,  there  will  be  one  and  only  one  solution  of 
equation  (3*7)  each  vcdue  of  n  between  1  and  n  Inclxislve,  and 
there  will  be  n  roots  of  equation  (2.10)  which  lie  on  the  real 
axis  between  u/c  and  u/c  . 

O  X 

As  might  be  anticipated,  we  cannot  solve  (3.7)  siq^licltly  for 
^  and  camnot  obtain  an  explicit  formula  Ibr  k^.  Nevertheless,  we 
can  solve  for  B  (or  u),  and  It  Is  there  that  Its  utility  lies. 

The  two  equations, 

B  -  w/w  »  [|(n-l)n+  d]/co8  n 


(3.9) 
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and 

k  -  («»/Co)(l  +  (3.10) 

vhlch  are  obtainable  from  (3*7)  and  (3*6) »  respectively,  relate 
k  and  ui  for  a  given  n  parametrically  throu^  Thiu  the  line 
k^(ui)  In  the  k>ui  plane  may  be  obtained  In  principle  by  allowing  ^ 
to  run  from  0  to  tt/2  and  computing  (u  via  (3*9)  and.  then  k  via  (3*10) 
for  each  valvie  of  u*  equations  (3.9)  and  (3.10)  therefore 
afford  a  complete  description  of  the  relationship  between  k  and 
(u  for  the  n-th  nomal  mode.  Furthermore,  as  shall  be  shown  In 
Section  V,  they  furnish  a  convenient  starting  point  for  analytic 
study  of  the  dispersion  characteristics  of  the  normal  modes. 

It  Is  to  be  noted  that  we  could  have  obtained  similar  equations 
even  if  we  had  not  assvmed  the  density  of  the  center  layer  to  be  the 
same  as  that  of  the  two  half  spaces.  With  the  same  deflation  of 
^  as  given  previously,  k  would  also  be  given  by  (3.10).  Instead  of 
(3.9)1  however,  we  would  find 

(n-l)TT/2  +  tan"^[(pypj^)tan  ^i] 

B  ■  '  '  » 

cos  u 

where  p^  is  the  density  In  the  half  spaces  and  p^  Is  the  density 

In  the  center  layer.  Ihe  relative  simplicity  achieved  by  equating 

p  and  p.  Is  clesu:. 

>'o  ^1 

Ihe  Complex  Roots 

The  study  of  the  roots  off  the  real  axis  may  be  greatly 
facilitated  by  use  of  the  auxiliary  equation  (3.9) »  although  the 
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resultlng  mathematical  expressions  are  much  more  complicated  than 
In  the  case  of  the  real  roots.  Here  ve  summarize  some  of  the  results 
of  such  a  study.  Die  mathematical  methods  are  discussed  In  Appendix 

A. 

The  complex  roots  In  the  first  q.uadrant  lie  along  a  line.  Die 
path  of  this  line  may  be  described  parametrically  by  eo^resslng  both 
the  real  and  Imaginary  parts  of  k  »  -t  1  k^  In  terms  of  the 
Imaginary  part  of  u*  This  line  of  roots  may  be  given  symbolically 
by  writing 

kj^  -  (m/c^)  f(k^cyu),  B,  A),  (3.11) 


since  a  unique  value  of  k^  may  be  obtained  for  every  value  of  k^ 
between  0  wd  <».  As  Indicated  in  the  above  expression,  the  path  of 
the  line  depends  on  B  and  A.  When  k^  Is  zero,  k^  will  be  greater 
than  u/c^.  For  sufficiently  large  B,  k^^  will  initially  decrease 
with  Increasing  k^  and  will  reach  a  minimum  which  is  between  0 
and  u/c^.  In  the  asymptotic  limit  of  large  k^,  the  path  Is  given 
by  the  relation 


-  (V%)A'V^log  [2(kjCyu))A"^],  (3.12) 


which  describes  a  slowly  increasing  function.  For  smaller  vedues  of 
B,  the  line  will  dip  only  silently  to  the  left  or  not  at  all,  and 
k^  will  always  be  greater  than  u/c^. 

Since  the  only  admissible  roots  are  those  lying  to  the  left  of 
the  branch  line  rising  from  '^ke  roots  of  (2.10)  In  the 
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flrst  quadrant  vlll  lie  along  the  portion  of  the  line  for  which 
kp  <  Obis  Implies  that  there  are  no  cooqplex  roots  for  small 

values  of  B,  and  that,  for  larger  values,  the  roots  all  lie  within 
a  finite  distance  of  the  real  axis— the  maxlmvm  value  of  being 
given  approximately  by 

kj  -  I  A  (3-13) 

We  Include  a  plot  of  the  line  of  roots  for  several  values  of 

2 

£  and  the  particular  case  A  «  0.2^  to  Illustrate  these  remarks 

u)  vs.  k^c^u)  as  we  wish  to  show 
the  path  of  the  line  relative  to  the  branch  line. 

As  In  the  case  of  the  real  roots,  each  complex  root  may  be 
associated  with  a  mode  number  m.  Ihe  m-th  root  will  lie  on  the 
intersection  of  the  line  of  roots  with  a  second  line, 

Sm^Vc/“''  ^3.14) 


(Fig.  4).  We  have  plotted  kj^c^ 


where  does  not  depend  on  B.  Illils  line  intercepts  the  real  axis 
at  some  point  between  uu/c^  and  o/c^.  For  large  values  of  k^  it  has 
the  asymptotic  form 


4 


(u)/Cj^)(mTT/2)(kj^cy-) 

log[2A"^pCyu)] 


(3.15) 


which  describes  a  line  of  slowly  decreasing  slope— k^  Increasing 
nearly  linearly  with  Increasing  k^.  Its  bdiavlor  for  Inteimedlate 
values  of  k^  Is  difficult  to  describe  by  analytic  methods.  How¬ 
ever,  the  path  of  the  line  may  be  found  for  any  desired  m  and  A  by 
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a  atzalc^tforvard  a\aierieal  conqputatlon.  In  Fig.  5  ve  show  the 

2 

results  of  one  such  calculation  for  several  values  of  m  and  A  «  .2^. 
It  Is  seen  that,  although  the  line  starts  and  ends  at  the  right  of 
the  branch  line,  It  dips  to  the  left,  so  that  Is  less  than  w/c^ 
for  a  certain  range  of 

Ihe  intersection  of  the  lines  (3*11)  and  (3.1^)  will  be  at 
most  in  only  one  point.  Ihey  will  intersect  only  if  B  Is  less  than 
some  value  vhlch  depends  on  m.  IThls  value  is  approximately  given 
by  (iih-1)tt/2,  vdilch  cuf fords  a  good  estimate  of  the  hlg^  frequency 
cutoff  for  the  coo^lex  modes. 

Ihe  intersection  of  the  two  lines  will  not  lie  to  the  left  of 


the  branch  line  xmless  B  is  between  two  limiting  values.  Biese 
two  values  may  be  obtained  In  principle  by  setting  k^^  equal  to 
u^c^  in  equations  (3«U)  and  (3>13)  an<i  then  eliminating  k^.  Tvo 
such  roots  for  B  will  always  exist,  provided  the  line  (3<^3)  Has 
in  part  to  the  left  of  the  branch  line.  These  two  values  of  B  may 
be  identified  as  the  high  frequency  and  low  frequency  cutoffs  (in 
units  of  w)  for  the  ra-th  mode. 

Another  nvnerlcal  scheme  may  be  devised  for  obtaining  the  varia¬ 
tion  of  k^2  and  kj^2  with  respect  to  B  for  each  mode.  Tbls 
scheme  Is  also  based  on  the  use  of  parametric  equations.  In 

o 

Figures  6  and  7  we  show  the  results  of  such  a  computation  for  A  ■ 
0.25  several  values  of  m.  The  low  and  high  frequency  cutoffs 
are  as  indicated. 


Our  study  thus  Indicates  the  following: 

1.  For  any  given  frequency  there  will  be  a  finite  nimiber  of 
roots. 
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2.  Ihe  roots  will  lie  along  a  line,  such  that  the  root  with 
lowest  mode  number  will  be  nearest  the  real  eucis. 

3>  With  Increasing  freq.uency  the  roots  with  low  mode  nunbers 
will  successively  disappear  and  roots  with  high  mode  nisnber  will 
successively  appear. 

4.  For  a  root  associated  with  any  given  m,  there  is  a  nonzero 
lover  bound  for  the  distance  of  this  root  from  the  real  axis. 
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IV.  THE  ASYMPTOTIC  SOLUTION 

Let  us  now  consider  the  nature  of  the  formal  solution  In  the 
event  that  r,  the  horizontal  distance  from  the  source,  is  large. 

Our  first  task  will  be  to  show  that  the  normal  node  solution  is 
the  dominant  portion  in  this  limit,  nils  is  not  immediately  clear 
for  our  particular  choice  of  source  and  receiver  positions,  EUid  seme 
consideration  should  be  given  before  concentrating  on  the  normal 
mode  solution. 

The  Direct  and  Reflected  Waves 

For  large  r,  the  net  distance  traveled  by  either  the  direct 
wave  or  the  wave  reflected  from  the  near  edge  of  the  center  layer 
is  nearly  the  same.  One  finds,  in  particular  that 

Rj  -  Rjj  -  2r"^(l2|  -  h/2)(|z^|  -  H/2),  (4.1) 

which  decreases  with  r  as  r  It  follows  that  the  direct  wave 

and  the  reflected  wave  will  arrive  at  virtually  the  same  time. 

Hie  time  lag  6t  of  the  reflected  wave  is  given  by  dividing  (4.1) 

by  c  .  Furthermore,  the  two  waves  will  have  almost  the  same 
o 

amplitudes,  that  of  the  reflected  wave  being  slightly  less  than 
that  of  the  incident  wave.  Since  the  two  waves  are  iSO*^  out  of 
phase,  they  will  very  nearly  cancel. 

The  sua  of  the  direct  and  reflected  waves  is  given  approximate¬ 
ly  by  the  ei^resslon 

t  Cf(T)  -  f(T-6t)]  +^f(T)  ,  (4.2) 


I 
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where  we  have  abbreviated  t  *=  t  -  Ihe  second  tern  de- 

.3 

creases  with  distance  as  r  and  Is  generally  negligible  at  large 
r.  In  the  event  the  function  f(T)  Is  differentiable,  the  first 
term  Is 

f(T), 

_2 

which  decreases  as  r  .  In  any  event  It  Is  clear  that,  for  any 
realistic  source  fxinctlon,  the  portion  of  the  wave  carried  by  the 
direct  and  reflected  waves  will  be  negligible  compared  to  any  wave 
which  decreases  as  r  ^  or  slower  at  sufficiently  large  r. 

The  Lateral  Wave 

The  lateral  wave,  which  Is  given  by  (2.14),  may  be  analyzed 
fl-4^ 

by  standard  methods.'  '  We  have  followed  the  method  attributed 
to  Gazaryan  which  Is  discussed  In  detail  by  Brehhovsklkh.^^^ 

It  Is  to  be  noted  that.  In  general,  U(B^)  Is  finite  at 
k  »  u/c^  and  hence  that  the  integrand  in  (2.l4)  vanishes  at  the 
beginning  of  the  path  for  the  k  Integration.  This  property  was 
achieved  In  the  early  stages  of  our  formulation  \dien  we  extracted 
the  direct  wave  and  the  wave  reflected  from  the  near  edge  of  the 
center  layer.  Since  the  Hankel  function  decreases  as  exp  (-k^r) 
with  Increasing  k^,  the  principal  contribution  to  the  integral 
will  come  In  the  earliest  portion  of  the  path  for  large  r.  If  one 
then  makes  approximations  consistent  with  this  observation,  he  may 
obtain  am  expression  \diich  is  proportional  to  r  . 

nils  procedure  Is  invalid  In  general,  however,  since  the  co¬ 
efficient  obtained  will  be  infinite  for  certain  discrete  frequencies 
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vhlch  correspond  to  the  cutoff  frequencies  of  the  noznsl  modes. 

When  one  considers  the  integration  for  such  a  frequency  at  the  out¬ 
set,  he  discovers  that  the  resulting  e:q>ression  is  proportional  to 

-1 
r  . 

A  general  expression  for  the  k  integration  vhlch  is  sufflcloatly 
valid  for  all  frequencies  may  be  derived  by  first  approximating 
U(0q)  “  U(“Pq)  points  on  the  branch  line  near  the  point  k  ■  od/c^ 
by  an  ei^ression  of  the  form 

11  ^ 

u(s„)  -  u(-s^)  . 


(where  T1  and  Q  are  independent  of  k)  and  then  integrating  over  the 
path  using  this  approximation.  When  one  does  this  and  assumes  r 
to  be  large,  he  may  obtain  the  following  expression  for  the  lateral 
wave: 


Re 


f «(») 


itt)(r/c  -t) 

i  (DB  sin  2B  -  cos  2B) 

vr 


dm,  (*^.3) 


where  E  is  the  frequency  in  units  of  i,  euad  D  is  defined  by  {2.5)* 
The  function  k(w2)  is  that  defined  by  Brekhovskikh  as 

4iW^  “  ”  "" 


k(w2)  *  itiit  r 

^/rT 


*  S^e”®  dS 


V7w"  ' 


ih.k) 


Its  argvment  is  dependent  on  r  and  m  and  may  be  taken  approximately 
as 

W^  »  I  r  AH"^  B  sin^  2B,  (4.5) 


I 
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vhere  A  is  defined  by  (3«3)>  derivation  of  (4.3)  le  lengthy 
but  straightforward  If  one  follows  the  st^s  outlined  by 
Brekhovsklkh. 

As  may  be  seen  from  the  asynptotlc  formulas  given  by 
Brekhovsklkh,  the  q.uantlty  K(W^)/W^  Is  of  negligible  magnitude 
when  is  large.  We  note  from  (4.3)  that,  for  large  r,  this  will 
be  true  unless  sin  2B  is  very  small.  Thus  the  significant  contri¬ 
butions  in  the  integration  over  tu  come  from  those  freq.uencies  near 
the  cutoff  frequencies. 

The  preceding  remarks  suggest  an  approximate  method  of  evaluat¬ 
ing  (4.3).  One  first  decomposes  the  integration  to  a  sum  of 
Integrals,  the  first  proceeding  from  ud  ^  0  to  in/4,  the  second 
from  ojn/4  to  mSrr/k,  etc.  In  each  of  the  terms  one  may  then  trans¬ 
form  the  variable  of  integration  to  W.  Since  the  range  of  W  will 
be  large  in  any  one  integral  and  the  principal  contribution  comes 
from  relatively  small  values  of  W,  the  limits  of  integration  in 
each  term  may  be  effectively  taken  to  be  0  and  <»  for  n  >  1,  and 
-®  and  ®  for  n  >  1.  Also,  since  the  factors  in  the  integrand 
other  than  K(W^)/W^  and  will  vary  slowly  with  W,  they 

may  be  replaced  by  their  values  at  W  >  0.  In  this  manner  one  may 
obtain  the  following  esgpresslon  for  the  lateral  wave: 

2  {si 


Z  -  „  1/2  - iw(r/o^-t)  . 

1  t(n-l)n/2]'^/2g(^^)e  o  F(T^)}  , 


tv>2 


(4.6) 
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vhere  m  ■  (n-l)Tiw  and 
c 


'it'')  ■  f  Att  ^ 

O  ^  * 


F(t)  »  J  x"^  K(x^)dx, 


All  the  Integrals  given  above  are  absolutely  convergent  and 
bounded  for  all  values  of  Tj^  and  T.  Furthermore,  for  any 
realistic  source  function,  f(t},  the  sum  In  (4.6)  is  absolutely 
convergent. 

Our  analysis  thus  shovrs  that  the  lateral  vave  Is  composed 

-4/3 

of  a  term  which  decree;  ies  as  r  '  and  a  term  which  decreases  as 

It  follows  that  the  lateral  wave  will  be  negligible  at  large 
r  when  compared  to  a  term  which  decreases  as  r"^. 

The  Complex  Mode  Solution 

The  conplex  mode  solution  is  composed  of  a  sun  of  Integrals 
over  u).  The  integral  for  the  m-th  mode  when  r  is  large  may  be 


written  in  the  form 


1  r*® 

hi 


(4.7) 


\diere  P  Is  a  bounded  function  which  does  not  depend  on  r.  Since 
m 
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the  Imaginary  part  of  k  increases  with  decreasing  frequency,  the 
major  contribution  to  the  Integral  comes  from  frequencies  close 
to  the  high  frequency  cutoff 

We  noted  in  the  last  section  that  generally  at  the  hlg^  fre¬ 
quency  cutoff,  the  imaginary  portion  of  k  Is  not  zero.  nils  Is 
true,  for  example,  in  the  cases  plotted  in  Fig.  7*  mlnlmun 
value  of  k^  vas  nearly  Independent  of  mode  number  with  a  value 
of  about  0.13(2/h).  nils  indicates  that  the  Integral  (4.7)  must 
decrease  eiqionentlally  with  r  and  hence  will  be  negligible  for 
large  r,  when  compared  to  any  term  decreasing  as  r”^  or  slower. 

It  follows  that  the  same  will  be  true  for  the  complex  mode  solution 
as  a  vdiole. 

nie  Normal  Mode  Solution 

nie  portions  of  the  wave  discussed  previously  have  all  fallen 
off  with  distance  faster  than  r”^.  As  is  well-known,  the  amplitude 
of  the  normal  mode  waves  fall  off  as  r'^  or  slower  (for  the  Airy 
phase),  nie  normal  mode  solution  is  therefore  the  dominant  portion 
of  the  solution. 

Evaluation  of  the  normal  mode  solution  depends  on  carrying 

out  the  Integrations  over  w  in  (2.13).  nils  Integration  may  be 

performed  to  a  good  approximation  by  the  method  of  stationary 
(l-4  2l) 

phase.'  ’  '  nils  method  rests  on  the  fact  that  the  dcmlnant 

contribution  to  the  Integral  at  large  r  comes  from  frequencies  near 
those  for  which 


Ijj  (kr  -  urt)  -  0 


(4.8) 
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or 

t  -  r/Vg,  (4.9) 

vhere  v  is  the  gro^p  velocity  defined  as 
S 

V  *  (^*10) 

o 

If  there  are  no  freq.uencle8  satisfying  (4.8),  the  integral  will 
generally  be  of  negligible  value  and  may  be  approximately  set  to 
zero. 

The  groi:^  velocity  depends  on  tu  and  the  mode  number  n: 

Vg  ■=  (^•^) 

In  the  next  section  we  shall  show  that  v  »  c  at  the  cutoff  fre- 

n  o 

quency,  that  it  initially  decreases  to  a  minimum  value  less  than 

c^  8uid  then  Increases  monotonically,  asymptotically  approaching 

Cj^  as  u)  -  ®.  This  is  identical  behavior  to  that  found  by  Pekeris^^^ 

for  his  model,  and  we  shall  follow  his  general  metnod.  We  shall 

denote  the  minimum  value  of  the  group  velocity  for  the  n-th  mode 

as  V  and  the  corresponding  V£dues  of  w  and  k  by  ui  and  k 
n,o  n,  0 

Similarly,  k'"  will  represent  the  third  derivative  of  k  (m)  with 
n,o  n 

respect  to  ui  at  ud  «  In  general  we  shall  use  primes  to  de¬ 

note  the  derivatives  of  k^^  with  respect  to  m. 

The  solutions  of  (4.8)  shall  be  denoted  u)  (t)  if  id  <  uu 

'  '  n,g'  '  n,o 

and  ID  (t)  if  w  >  (u  .  where  t  <■  c  t/r.  (Die  subscripts  g  and 
n,\r'  '  n,o'  o''  o 

w  stand  for  "ground"  and  "water,"  respectively,  and  were  chosen 
in  order  to  make  use  of  the  temlnology  Introduced  by  Pekerls.) 
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Recalllng  our  previous  remarks  about  the  behavior  of  the  groxqp 

velocity j  ve  see  that  u)  Is  defined  for  all  t  between  1  and 

n,g 

T  -  ■  c  /v  while  u)  Is  defined  for  all  t  between  c/o, 
n,f  o'  n,o'  n,w  o'  1 

and  T 

n,f 

By  use  of  the  method  of  stationary  phase^  a  typical  Integral^ 

'f’n  *  J  ^  e“^“^  duj,  (^-12) 

appearing  In  (2.13)  becomes,  for  sufficiently  large  r; 


0 

T  <  1; 

(^.13a) 

'*’n,g' 

1  <  T  <  cjc^i 

(^.13b) 

'*’n,g  *  ^n,v* 

c 

0 

—  <t<t  --e; 

Cj^  n,f  n' 

(4.13c) 

\,Ai* 

T  >  T  _  -  G  ; 

n,f  n' 

(4.13d) 

where  the  quantities  cp^  ,  cp  and  'o  are  generlcally  called 

n,g'  ^n,w  n,Al 

the  "ground  wave, "  the  'Sraiter  wave, "  and  the  "Airy  phase, " 

respectively.  Bie  choice  of  the  parameter  is  somewhat  arbitrary — 

the  prime  restriction  Is  that  cp  „  9  approximately  equal 

n^v 

to  fp^  when  ’’’  *  ‘’’n  f  "  ^  sufficiently  large.  A  satis¬ 

factory  choice  Is 

«„  ■  “  '0(5 


The  mathematical  expressions  for  9  and  9  are  as 

ll^W 


follows 
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,  ,  i(k  r-uut), 

fn.g  '  r  {'  }  , 

(U  ■=  U)  (t) 

a,g'  ' 


\rtiere 


“■  «n,g<’) 


U)  =  U)  (t) 


cmd 


.  ,  i(k  r-urt). 

'Pn^v  =  7  {«  } 


(^•15) 


(1^.16) 


(4.17) 


U)  =  U)  (t) 

n,g^  ‘ 


(4.18) 


The  quantities  and  are  defined  analogously  to  (4.l6) 
and  (4.17).  Wie  quantity  «  . .  is  given  by  the  emression 

n^Ai 


^p, 


iTT/4 

'n,Ai  “  y5/b 


g(tt)  )R(u> 
n,o'  '  n,o 


A^(X),  (4.19) 


n.o 


where 


R(u;  ) 
'  n,o' 


(4.20) 


and  Al(X)  is  the  Airy  function  defined  as 

Ai(x)  -  (2n)"^  J  *“^du.  (4.21) 


Its  argument,  X,  is  given  by 
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X  =  [t  -  (r/v^^^)]  (|k”;Jr/2) 


-1/3 


(4.22) 


nie  method  ve  have  Indicated  for  evaluating  the  Integral 

(4.12)  1&  approximate  and  Is  generally  regarded  as  Inadequate 

for  times  near  r/c^  and  r/c^.  However,  for  our  particular  model, 

cp„  vanishes  at  t  =  r/c  and  m  vanishes  at  t,  ■  r/c, .  We 
^n,g  '  o  ’^n,w  1  '  1 

accordingly  feel  that  the  solution  (4.13)  gives  a  qualitatively 

correct  value  of  the  Integral  for  the  times  when  Its  value  Is 

largest.  It  will  not  give  a  true  value  of  the  amplitude  of  the 

first  peak  in  the  ground  wave,  but  this  will  be  of  negligible 
magnitude  at  large  r,  even  If  one  uses  a  more  refined  method. 

With  this  approximation  the  normal  mode  solution  becomes 


♦ 


I»1 


n=l 


Re  tpjj(r,t). 


(4.23) 


Itie  Interpretation  of  this  solution  will  be  our  principal  goal 
In  the  remainder  of  this  paper. 
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DISPBRSIVE  CHARACTHUCBTICS  OF  THE  NORMAL  MODES 


nie  dispersive  characteristics  of  the  n-th  normal  mode  are 
those  quantities  and  functions  derivable  from  which  are 

needed  for  the  evaluation  and  Interpretation  of  the  quantity  (p^Cr^t) 
as  given  by  eqmtlon  (4.3)  In  the  preceding  section.  A  partial 
list  of  such  dispersive  characteristics  may  be  given  as  follows: 

1.  nie  phase  velocity  v^(ui)  »  lu/h 

2.  The  group  velocity  v^(u)) 

3.  The  minimum  value  v  of  the  group  velocity 

n^o  -TT 

4.  The  frequency  u)^  ^  at  which  the  group  velocity  is  a 
minlmvBi 

5.  The  second  derivative  h^(u>)  of  k  with  respect  to  u) 

6.  The  third  derivative  k"'  at  u)  »  u)„  „ 

n^o  n^o 

7.  The  functions  m  (t)  and  u)  (t) 

n,g  n^w 

8.  The  quantity  t  _  which  represents  the  time  of  arrival 

n,r 

of  the  Airy  phase  in  units  of  r/c^. 

These  dispersive  characteristics  are  Inherent  properties  of  the 
model  and  are  Independent  of  the  positions  of  the  source  and  re¬ 
ceiver.  In  this  section  we  shall  study  these  chsuracterlstlcs  in 
some  detail— with  particular  emphasis  on  their  depmidence  on  mode 


nvnber. 


The  function  k^(w)  is  given  by  the  parametric  equations  (3.9) 
and  (3.10).  As  ^  increases  from  0  to  n/2  both  k  and  lu  increase 
monotonlcsdly.  It  follows  that  k^(<u)  is  a  monotonlc  increasing 
function  of  w. 
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It  Is  possible  In  principle  to  obtain  by  first  solving 

(3*9)  for  u  In  terms  of  w  and  then  siibstltutlng  this  Into  (3.IO). 
Sils  Is  rigorously  valid  as  there  is  a  one  to  one  relationship 
between  g,  and  w  over  the  range  of  interest.  We  may  therefore  write 
k^(uu)  formally  as 

ItnC")  *  T(g^(B)  ) ,  (5 .  l) 

provided  It  Is  remembered  that  g  euid  are  defined  only  for 
B  >  i(n-l)TT.  Here  we  have  abbreviated 

T(n)  »  (1  +.  Bln^u)^/^. 


The  phase  and  group  velocities  are  then  given  formally  by 
the  expressions, 

Vp  =  c^T  (5.2) 

and 

Vn  -  [T  +  B(dT/dti)(d4^/dB)]‘^,  (5.3) 


\diere  g  Is  to  be  evaluated  at  g  =  derivative  dg^/dB 

In  (^.3)  may  be  e]q;)reBBed  in  terms  of  g  since 


dg  ^dB\-l  cos  g 
dB  “  'dg'  “  1  -I-  B  sin  g  ' 


(5.4) 


by  differentiation  of  (3.9)*  With  the  use  of  this  expression,  we 
nay  obtain  the  following  expression  for  the  group  velocity: 
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V 

n 


B  +  Bln  u 
1  +  B  sin 


1  + 


.2  , 

A  sin  p, 


(5.5) 


The  group  velocity  corresponding  to  any  desired  frequency  above 
cutoff  may  be  obtained  by  evaluating  (5*5)  at  p,  =  P^(b). 

An  alternative  point  of  vlev  is  to  regard  (3.9)  and  (^.2)  as 
parametric  equations  for  the  i^se  velocity,  and  (3.9)  and  (^.^) 
as  parametric  equations  for  the  group  velocity.  Bie  functions 
Vp(uj)  and  v^(<*')  ®ay  be  tabulated  by  letting  p  range  from  0  to 
tt/2  and  computing  lu  =  B  i  from  (3*9)  and  then  v^  and  v^  from 
(5.2)  and  (5.5)  for  each  value  of  p. 

In  a  manner  similar  to  that  used  for  deriving  parametric 

eqviatlons  for  v^(<“)#  we  may  derive  parametric  equations  for  the 

derivatives  (of  any  order)  of  k  or  v  with  respect  to  frequency. 

n  n 

In  each  such  case,  we  may  express  the  desired  dependent  variable 
as  a  function  of  p  and  B,  which  may  be  written  e:q)licitly  in 
terms  of  elementary  functions.  Such  em  equation,  in  conjunction 
with  the  equation  (3.9)  for  u)  in  teims  of  p,  gives  a  parametric 
representation  of  the  variation  of  the  dependent  variable  with  ui. 

One  such  expression  which  will  be  useful  in  our  analysis  is 
that  giving  the  first  derivative  of  the  group  velocity.  We  find 


dv  (u() 
du 


-  [(x-B)®  -  f?  ■  1|],  (5.6) 


n 


Sr^B^c 
p  o 


B^  •  (l+B  sin  p)ycos  p. 


where 


(5.7) 
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X  «  B  COB  M.  “  ^(n-l)TT  +  n,  (5*8) 

0  =  ^(l+A^)  ^  cos  |ji  Bin”^  li  (l+A^  BinV)>  (5*9) 

T1  =  2(1+ A^)"^  COB^^,  (1+A^  Bin^p,),  (5*10) 


and  is  given  by  (5*5) •  Expressions  for  and  k^  may  be  found 
directly  from  (5-5)  and  (5 •6),  since 


and 


k* 


n 


k* 

n 


-(v  ) 
'  n' 


-2 


(dv^/du)). 


(5.11) 


(5.12) 


Ei^jressions  for  higher  order  derivatives  become  more  complicated. 

We  have  carried  this  as  far  as  the  seventh  order  but  will  not 

reproduce  the  results  here.  The  only  higher  order  derivative 

needed  is  k'"  ,  whose  value  at  uu  determines  the  amplitude  of 
n  n,o 

the  Airy  phase.  Our  expression  for  this  derivative  is  given  in 
Appendix  B. 

The  Phase  Velocity 

The  phase  velocity  as  given  by  (5.2)  is  a  monotonic  decreas¬ 
ing  function  of  pi.  Since  pi  increases  with  frequency,  the  phase 
velocity  will  also  decrease  with  iw.  At  pi  =  0  the  phase  velocity 

is  c  ,  while  it  is  c,  at  pi  -  It  follows  that  v^  decreases 
o'  12  p 

monotonically  from  c^  to  c^^  as  the  frequency  increases  from 
to  ».  For  frequencies  very  close  to  the  cutoff  frequency  one  may 
derive  a  power  series  expansion  by  finding  the  derivatives  of  v 
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vith  respect  to  u)  In  texns  of  )ji  and  then  setting  ■  0.  In  this 
manner  one  finds 

V®o  •  ^  ‘  I  +  •••>  (5-13) 

where  ■  i(n-l)n  and  6B  »  B-B^.  In  the  other  limit  of  very 
large  frequencies  one  may  obtain  an  asymptotic  expression  for 
by  expanding  It  In  terms  of  B~^.  This  Is  accomplished  by  expending 
B  ^  In  terms  of  (n/2)-^  using  (3*9) >  and  then  inverting  this  series 
and  substituting  the  resulting  esgpresslon  Into  (3*2} •  Doing  this, 
ve  find 

Vp/cj^  =  1-  I  A^(A^+l)"^(|nn)V^ 

By  using  these  two  limiting  expressions  and  because  v^  decreases 

with  frequency,  one  may  obtain  a  fairly  accurate  sketch  of  the 

function  V  (w). 

P'  ' 

The  Group  Velocity 

Analogous  techniques  may  be  applied  to  the  study  of  the  grovp 
velocity  For  frequencies  near  cutoff  it  behaves  as 

V<=o  “  ^  *  i  A^(2b2a2+3b2.3)(6b)2,  (5.15a) 

while  for  large  frequencies  it  behaves  as 

V°1  •  ^  I  A^(l+A^)“^(|nn)V2. 


(5.15b) 
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It  iB  clear  from  these  expressloiiB  that  at  the  cutoff 

frequency  and  that  it  initially  decreases  with  Increasing  frequency. 
For  large  frequencies  is  less  than  c^  emd  Increases  with  increas¬ 
ing  frequency— approaching  c^^  as  ou  -•  <».  it  follows  that  must 
have  a  minimum  value  somewhere  which  is  less  than  c^. 

Biat  has  only  one  minimum  can  be  verified  by  studying 

equation  (5*6).  It  is  readily  seen  that  dv  /diu  is  positive  or 
negative  if  the  quantity  (x-a)  -  0  -  11  is  positive  or  negative 

since  all  the  other  factors  are  positive.  !Biis  quantity  in  turn 
will  be  positive  or  negative  if 

X  -  a  -  (a^+'n)^'^^  (5.16) 


is  positive  or  negative.  The  quantity  x  Increases  with  p,  with  a 

slope  of  unity.  The  quantity  a  +  (S  +7))  '  is  infinite  at  »  0 

and  will  generally  decrease  with  Increasing  its  veilue  being 

zero  at  u  =  For  large  A  it  may  have  a  positive  slope  over  a 

small  range  of  u*  maximum  positive  slope  it  can  have,  hc^wever, 

is  less  than  l/2.  (ihis  may  be  verified  by  a  careful  examination 

of  the  expressions  for  a  and  T).)  It  follows  that  the  quantity 

(^.l6)  is  a  monotonlc  increasing  function  of  )ji.  It  is  minus 

infinity  at  pi  -  0  and  m\/2  at  =  tt/2.  It  therefore  must  be  zero 

in  one  and  only  one  place  at  some  nonzero  value  u  .  It  follows 

n ,  0 

that  dv^/du)  is  less  than  or  greater  than  zero  depending  on  whether 

u)  is  less  than  or  greater  than  lu  where  (u  is  that  obtained 

n.o'  n.o 


from  (3-9)  with  j*  « 


1 


1 
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In  Fig.  6  ve  give  a  schematic  sketch  of  with  the  | 

characteristic  behavior  predicted  above.  The  actual  detailed  shape 
2 

Is  dependent  on  A  and  n.  A  good  quantitative  specification  of 

the  shape  of  the  group  velocity  curve  may  be  attained  by  focusing 

one's  attention  on  a  relatively  small  number  of  parameters  which 

characterize  the  curve.  These  are  (1)  6B  ,,  the  frequency  relative 

n,l 

to  cutoff  (in  units  of  ou)  at  idilch  the  group  velocity  Is  equal  to 
(2)  6B^  the  frequency  relative  to  cutoff  at  which  the  group 
velocity  Is  a  minimum,  (3)  v^  the  minimum  value  of  the  group 
velocity,  and  (4)  v^  the  curvature  at  the  grovq?  velocity 
minimum.  (The  first  three  of  these  parameters  are  Illustrated 
In  Fig.  8*)  These  four  parameters,  with  the  two  equations  (5*13) 
and  (5.14),  give  a  fairly  complete  specification  of  the  function 
v^(u)).  The  alternative  would  be  to  tabulate  v^(u))  numerically. 

This,  however,  would  be  a  cumbersome  way  of  depicting  the  varia¬ 
tion  of  the  groiqj  velocity  curves  with  mode  nuniber. 

One  finds,  by  equating  (5*5)  to  c, ,  that  the  parameter  6B  , 

X  11^  X 

corresponds  to  a  value  .  which  satisfies  the  equation 

n,i 

(|(n-l)Tr  +  u)tan  u  -  C(l+A^  sin^u)/(l+A^)]^^^-  (5.17) 

For  large  n,  the  solution  of  this  equation  Is  asymptotically 

6  -  (l+A^)’^/2(i(n-l)TT)"S  (5.18) 

the  solution  for  smaller  values  of  n  must  be  found  numerically. 

Since,  for  small  |ji, 

I 

I 
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2 

6B  ~  p,  +  B  p  fZ, 


we  find 


6B 


"  (|(n-l)TT)(l+A^) 


(5.19) 


for  large  n.  A  plot  of  6B  .  versus  n  is  given  in  Fig.  9  for  the 

n^x 

2 

particular  case  of  A  =  O.5625  (or  “  0.0).  The  shape  of  the 

curve  for  small  values  of  n  was  found  by  computing  v^(u))  and 

graphically  detenalning  6B  . . 

n,j. 

The  characteristics  of  the  groiq?  velocity  curve  at  its  minimum 

may  be  studied  analogously.  The  value  of  p  may  be  found  by 

n^  0 

solving  the  eq.uatlon 


i(n-l)TT  =  -  p  +  0  +  (B^+ 


(5.20) 


where  0  and  T)  are  given  by  (5.9)  and  (5.I0).  For  large  n,  one  finds 


=  r_  2  1^/3 

W-|..2\  /  ,\J  ' 


(l+A  )TT(n-l)‘ 


(5.21) 


and,  therefore. 


®®n,o  *  I  (l+A^)"^/^, 


(5.22) 


and 


r.o  c 


£(A^+l)®^3(i(n-l)TT)2^3 


). 


(5.23) 


1 


1 
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An  asyn^totlc  expreeslon  for  v”  ^  may  also  be  derived  by  first  ex¬ 
pressing  v"  in  tezms  of  ^  and  then  setting  u  -  m>  .  and  taking  the 
n  0 

limiting  fozu  for  large  n.  We  find 


n.o 


w  ^  (|(n-l)TT)^^^ 


(5.24) 


A  simple  scheme  may  be  used  for  finding  the  variation  of  the 

quantities  6B  .  v  and  v"  with  n  which  avoids  the  computation 
n^o  n,o  n>o 

of  the  complete  curve  for  each  value  of  n.  Cne  lets  u  run  over  a 

sequence  of  values  from  tt/2  to  0.  For  each  he  computes  n  from 

(5.20).  He  then  uses  this  n  to  compute  B  and  B  from  (3.9)  and 

c 

finds  the  difference  6B  ..  Then  he  uses  these  values  of  B  and  u 

n>o 

to  find  v^  and  v".  The  quantities  6B,  v^,  and  v^  when  plotted 
versus  n  give  the  behavior  of  the  desired  characteristics  with 
mode  nunber. 

In  Figures  10,  11  and  12  we  show  the  results  of  such  a  calcula- 
2 

tion  for  A  ■  O.5625.  Figure  10  shows  the  minimum  groqp  velocity 

decreasing  monotonlcally  with  increasing  n,  approaching  c^/c^  as 

n  approaches  infinity.  Hie  quantities  6B  and  v"  do  not  behave 

n,o  n,o 

as  systematically.  We  see  that  6B  has  a  mlnlmvai  value  at  n  ~  6 

Qy  O 

and  that  v"  has  a  maximun  value  at  n  — 
n,o 

The  quantity  k'*'  may  be  determined  from  v"  ,  since 
n,o  n,o' 


k. 


n,o 


-4/3 


For  large  n  this  is  proportional  to  n 
(5.23)  and  (5.24). 


as  may  be  seen  from 
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Varlatlon  of  k"  With  Freguency 

The  analyela  of  the  q,uantlty  may  be  accomplished  by  methods 
similar  to  those  used  in  studying  the  group  velocity.  Since  it  is 
related  to  the  derivative  of  the  group  velocity  by  the  relation 
($.12)4  it  is  dlear  that  it  must  have  one  of  the  two  geieral  fozms 
depicted  in  Fig.  I3— depending  on  whether  n  >  1  or  n  >  1. 

If  n  =  1)  is  zero  at  uu  »  0  and  initially  increases  with 
Increasing  frequency,  such  that 

k"  -  m  B  [3  -  (5/6)(3A^+16)B^]  (5.25) 

for  small  B.  It  reaches  a  maximum  value  \dilch  is  of  the  order  of 
magnitude  of 

2  w 

S  (3A  +16)^^^ 

and  then  decreases,  reaching  zero  at  (u  >  It  continues  to 

decrease  to  some  minimvm  value  auad  then  Increases  again,  its  asymptotic 
fozm  being 

kj|  —  -  Att)^/B^.  (5-26) 

If  n  >  1  the  quantity  k^  has  a  nonzero  positive  value  at  the 
cutoff  frequency  and  decreases  llneeurly  until  it  reaches  a  mlnlmun 
value  which  is  negative  and  then  increases  again,  its  asynqirtotlc 
fona  being 


(5.27) 


-U2. 


(|ArxTT)2/B3. 

We  shall  be  interested  primarily  in  the  behaviour  of  k^'  for 
frequencies  less  than  cu^  For  large  n  the  behaviour  over  this 
range  of  frequencies  is  given  approximately  by  the  formula, 

K  =  - - ^  (1  -  (1+A^)B  (5.28) 

c  (uCl+B  6B]-^/  ° 

o  c 

idiere  B^  is  l/2(n-l)n  and  &B  is  B-B^.  We  may  see  from  this  equation 
that  k^  is  large  when  6B  is  zero  but  decreases  rapidly  with  in¬ 
creasing  frequency.  If  6B  is  much  greater  than  B^^,  k^  is  decreas¬ 
ing  as  6B’^^^.  When  6B  is  zero,  k"  is  approximately  proportional 

-I/2 

to  n  and  when  6B  is  unity,  k^  is  proportional  to  n  '  .  Ihe  conse¬ 
quences  of  this  type  of  dependence  on  n  will  be  discussed  in  the 
next  section. 

Variation  of  the  Dispersion  Characteristics  With  t 

The  parameter  t  represents  time  in  units  of  r/c^  and  is  formally 
related  to  v^  by  the  equation 

T  -  (5.29) 

Let  us  note  that,  using  (5*5)/  ve  may  also  express  t  in  terms  of  p,* 
Thus  the  variation  of  any  characteristic  with  t  may  be  found  in  the 
manner  that  one  would  find  its  variation  with  frequency.  One  ex¬ 
presses  both  T  and  the  dependent  variable  in  terms  of  p,  and  has  two 
parametric  equations  which  may  be  used  for  either  nmerlcal  computa¬ 
tions  or  as  a  basis  for  deriving  approximate  ejqpresslons. 
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Ue  summarize  here  some  of  the  properties  of  the  two  functions 
(I)  and  (JO  .  Most  of  these  properties  may  he  readily  derived 

rijV 

from  the  preceding  discussions  In  this  section. 

For  T  between  1  and  c  /c, ,  ou  Increases  from  the  cutoff 

o'  n,g 

frequency  to  (jj  , .  For  t  close  to  1  we  find 
n,± 

“n^g  *  *c  (5*30) 


for  n  ■>  1,  and 


“n^g  ”  “  (2/3)^/^  a'^/^(t-1)^/^ 


(5.30b) 


for  n  =  1.  The  frequency  thus  varies  slowly  with  t  in  this  limit 
for  large  n. 

For  T  near  u;^  ^  varies  as 


(Co/ci)i  A(|njT) 


(5.31) 


If  T  is  close  to  T  we  have 
n,f' 


=  u)  +  f  ..'TiTs  " 


(5.32) 


where  t  _  is  c  /v  .  One  may  study  the  variation  of  this  draendence 
n  ^  r  o  n  ^  o 

with  n  by  using  the  asymptotic  formulae  (5*23)  and  (5*24). 

O 

The  quantity  ^  is  asymptotically  .  For  large  n  it 

*  -2/3 

differs  from  this  by  a  quantity  proportional  to  n  ^  ^  as  may  be 
seen  from  (5.23). 
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We  shall  also  need  to  kaov  hov  the  phase  velocity  of  the  ground 
wave  varies  with  t.  When  t  ■  1,  the  phase  velocity  Is  c^.  For 
o  ■  proportional  to  (t-i)  for  t  close  to  Ij  for  larger 

O 

n.  It  Is  proportional  to  (t-1)  .  The  following  eq^uatlon  gives  a 
reasonable  approximation  for  t  between  1  and  ^  if  n  is  large: 


V  »  c 
p  o 


[-1 


(t-1) 
B?(1+A^ 


(5.33) 


(The  derivation  Is  discussed  In  Appendix  C. )  It  Is  clear  that  v 

P 

decreases  relatively  slowly  with  t  If  n  (or  B^)  Is  large. 

One  may  readily  show  from  (5*33)  that,  \Aen  t  =  t  «,  the  phase 

n,r 

velocity  differs  from  c^  by  a  (juantlty  proportional  to  n"^' ^  in  the 
limit  of  large  n.  It  follows  that  the  phase  velocity  of  the  ground 
wave  is  nearly  constant  In  the  limit  of  large  n. 


VI.  WB  E3CCI1AT10W  AMPLETUPBS 

The  noxtnal  mode  solutlos  given  In  Section  IV  may  be  interpreted 
as  follovs:  The  waveform  at  time  t  ■  rr/e  consists  of  a  sum  of 
waves  of  slowly  varying  fre<iuency.  The  freqiiencies  present  are 
u)^  ^(t)  and  ^(t)  for  n  «  1^2, 3>  •••  •  In  the  previous  section, 

we  showed  how  one  may  determine  Just  what  frequencies  correspond  to 
a  given  t.  Otiese  frequencies  were  independent  of  the  location  of 
the  source  and  of  the  time  dependence  of  the  source. 

The  next  question  we  consider  is  that  of  the  relative  aiiq)lituides 
of  the  various  waves.  It  is  clear  from  equations  (4.1^  -  4.21) 
that  these  amplitudes  depend  on  a  nunber  of  factors.  One  of  these 
is  g((u),  the  freqviency  spectrum  of  the  source.  Here  we  shall  con> 
slder  the  factors  other  than  g(w),  l.e.,  those  conqprlslng  the 
excitation  amplitude  for  that  particular  frequency. 

nre  excitation  amplitude  for  a  particular  frequency  is  given 
by  different  expressions,  depending  on  whether  the  frequency  is 
associated  with  the  ground  or  water  waves  or  with  the  Airy  phase. 

We  define  the  following  quantity: 


ifu)  <  w  <  w  „-a„oru)>m^  +  cr^,  or 
c  n,o  n  n,o  n 

.  (0.536)r*^‘ 

(6.1b) 

if  I*  -  <  V 


where 

,  (6.1c) 

as  the  relative  excitation  amplitude  of  the  n-th  normeil  mode.  It 
is  a  dimensionless,  real,  and  positive  quantity  which  may  be  used 
to  determine  the  relative  strengths  of  excitation  for  various  fre¬ 
quencies.  In  conjunction  with  some  knowledge  of  the  source  spectrum, 
the  knowledge  of  the  relative  excitation  functions  for  each  normal 
mode  will  allow  us  to  make  some  qualitative  predictions  of  the 
waveform.  In  particular,  we  may  predict  Just  which  frequencies 
associated  with  \dilch  normal  modes  will  be  dominant  at  any  given 
time,  and  we  may  predict  the  times  at  \dilch  the  wave  amplitudes 
will  be  greatest. 

We  shall  not  concern  ourselves  with  the  relative  phases  of  the 
various  normal  modes.  Ihese  will  vary  erratically  with  distance. 

It  is  more  instructive  to  examine  those  features  of  the  wave  which 
vary  slowly  with  distance  at  large  distsuices.  For  qualitative 
estimates  we  have  replaced  the  Airy  function  by  its  maximum  value, 
.536,  in  (6.1b).  This  will  suffice  when  comparing  the  relative 
magnitude  of  the  Airy  phases  to  other  portions  of  the  wave. 

TOie  Fvmctions  G(w) 

All  dependence  on  the  source  and  receiver  positions  in  the 
relative  excitation  functions  is  contained  in  the  functions 
G(k^)  »  G^(u)),  which  we  defined  in  (2.12).  If  we  use  the  definitions 
of  u  in  equations  (3-4),  the  quantity  may  be  written  as 
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G^(u,) 


B  sin  u 
1  -«■  B  Bin  pi 


COB  pi  e 


-DB  Bln  p, 


(6.2) 


Since  B  may  also  be  expressed  In  texins  of  pi,  the  equations  (3*9) 

and  (6*2)  constitute  a  pair  of  parametric  equations,  and  we  may 

numerically  tabulate  versus  uu  In  the  manner  described 

previously  for  the  dispersion  characteristics. 

If  one  were  to  carry  out  such  a  calculation  he  would  find  that 

G^(u))  has  the  schematic  form  shown  In  Fig.  l4.  It  Is  zero  at  the 

cutoff  frequency,  rises  to  a  single  maximum  G  at  6B  >  d  and 

n,o  n,o 

then  decreases,  asymptotically  approaching  zero  with  Increasing 
freq\jency. 

a 

Analysis  of  the  parametric  equations  shows  that  G^  behaves  as 

G^((b)  *>  B^  n  ■  1  (6.3a) 

-  B^6B  n  >  1  (6.3b) 


for  frequencies  near  cutoff,  emd  as 


1  -2  ”(®  ”  o*^)® 

Gn('")  “  (2^^)  ®  ® 


(6.4) 


for  very  large  frequencies.  Die  behaviour  for  intermediate  fre¬ 
quencies  may  be  studied  by  concentrating  on  the  four  quantities 

G  ,  A  1  A  , ,  and  A  irtilch  give  the  maximum  value  of  G„, 
n,o'  n,o'  n,l*  n,2'  ®  n' 

and  the  frequencies  relative  to  cutoff  In  units  of  lu  at  which  G^ 

Is  maxlmua  or  l/2  Its  maximum  value.  (See  Fig.  l4.)  Die  difference 

A  --  A  .  ■  BW  defines  the  band  width  carried  by  the  n-th  noznal 
n,e  n,i  n  — 


mode. 
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For  large  n,  the  quantity  „  is  Independent  of  n  and  coin- 

O 

cldes  with  the  maximiai  of  x(l+x)”^e"^.  For  D  -  1,  this  maximum  is 
about  0.2.  For  larger  0  It  decreases  as  while  for  smaller  D 
It  increases j  and  asymptotically  approaches  1  for  very  small  D. 

In  this  limit  we  find  that  and  A^^  g  inversely 

proportional  to  B^(or  n);  the  coefficients  will  depend  on  D.  For 
D  ■  Ij  we  find  the  following: 

^n,o  * 

An^l  =  0.15  (6.5b) 

Va  =  3.53  b;\  (6.5c) 

«  3.38  B”^.  (6.5d) 


IBie  coefficients  will  decrease  roughly  as  D  ^  for  larger  values  of  D. 

Another  quantity  of  interest  is  the  value  of  G  (<t»)  vdien  u;  >  w 

n'  '  n,o 

If  we  substitute  (^>21)  suid  (5*22)  into  (6.2),  we  find  for  large  n 
that 

Gn(<«n,o^  =  C-  -  (l/6)(l+A^)’^],  (6.6) 


which  decreases  rapidly  with  increasing  n.  In  a  similar  manner,  we 

may  find  G_(u)  ,)  where  uu  ,  is  the  frequency  at  which  the  group 
n  n,l  n,i  ^ 

velocity  is  c^.  We  find,  in  the  limit  of  large  n,  the  following 
expression: 


In  this  limit,  Gq(<i^^j^)  Is  independent  of  n. 


(6.7) 


-1^- 


Ihe  Quantity  ^ 


Let  us  define 


E°  («.) 


c;\(u.) 


172 


(6.8) 


nils  Is  eq.ual  to  the  relative  excitation  function  everyvhere  except 


In  the  vicinity  of  (u^ 

nils  quantity  may  also  be  e^g^tressed  In  terms  of  ^  and  nay  be 
tabulated  by  letting  through  a  sequence  of  values  and  comput¬ 

ing  (u  and  Cr^(u))  for  each  llhe  schematic  form  of  this  function  Is 
Illustrated  In  Fig.  1^.  It  Is  zero  at  the  cutoff  frequency.  For 
larger  frequencies  It  Increases  to  a  maxlmvan  and  then  dec]*eases. 
Bids  decrease  is  interrupted  by  a  narrow  singularity  at  For 

frequencies  greater  than  tu  .  £°  decreases  rapidly,  asymptotically 

Q 

approaching  zero. 

One  may  study  the  function  E°  by  methods  similar  to  those  used 
previously.  Our  analysis  reveals  the  following: 

1.  For  very  large  frequencies,  the  asymptotic  form  Is 

-  D(B-  Inn) 


E°  «=  ^  e 
n  AB 


(6.9) 


2.  In  the  limit  of  large  n,  the  magnitude  of  E°  at  Its  first 

maximum  decreases  with  n  as  B  For  D  >  1,  the  coefficient  Is 

c 

approximately  0.^3  A 

3.  Bils  maxlmvan  occurs  at  a  value  of  6B  \dilch  in  the  limit  of 
large  n  Is  proportional  to  b”^. 
limiting  form  to  be 


For  the  case  D  «  1,  we  find  the 


1 


I 
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6B  >  2.10  B 


-1 


(6.10) 


4.  Ibe  mlnimun  value  of  B  between  Its  maxtmm  and  the 

n 

singularity  at  lu^  ^  Is  very  small  and  decreases  maxicedly  with  In¬ 
creasing  mode  number.  In  the  limit  of  large  n,  we  find  the  value 
of  the  function  at  this  minimum  to  be 


aCi+a*^) 


13^^  exp 


5.  The  minimum  occurs  when  lu  is  very  close  to  tu^  For  very 


large  n,  we  find  the  difference  to  be  proportional  to  B 


•V3. 


u> 


n.o 


2B^/3jj 


2sV6 


(6.12) 


Amplitude  of  the  Airy  Phase 

Ihe  relative  excitation  function  £^(tt))  for  the  Airy  phase 
(u)  -  u)^  ^)  Is  proportional  to  r^^  and  Increases  slowly  with  r. 

It  Is  clear  that,  for  very  large  r,  the  amplitude  of  the  Airy  phase 
will  be  greater  than  that  of  the  maximum  amplitude  of  the  ground 
wave.  It  Is  Instructive  to  see  Just  how  large  r  must  be  In  order 
that  this  be  true. 

We  set  J^(Wj  q)  «l«al  to  evaluated  at  the  frequency  for  which 
the  ground  wave  amplitude  Is  the  largest.  This  gives,  for  D  ■  1, 

In  the  limit  of  large  n,  the  following  approximate  value  for  r: 


r 

H 


A 


1 
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To  get  an  Idea  of  the  magnitude  of  r,  ve  consider  the  case 
n  «  10,  A  •  0.25*  This  gives  r/H  of  the  order  of  10  .  We  can 
conceive  of  no  practical  case  in  which  the  ratio  would  ever  be  this 
large.  If,  for  example,  one  were  to  try  to  apply  this  model  to 
sound  propagation  In  the  atmosphere,  H  would  have  to  be  of  the  order 
of  20  km,  and  r  would  be  of  the  order  of  500  times  the  circumference 
of  the  earth. 

^Is  demonstrates  that,  unless  one  Is  Interested  in  very  low 
freqviency  sources,  the  dcmilnant  waves  carried  will  be  the  ground 
waves.  We  shall  henceforth  concentrate  our  attention  on  these  waves. 
Time  Dependence  of  Wave  Amplitudes 

Having  eliminated  the  Airy  phases  and  water  waves  from  our 
consideration,  we  shadl  consider  only  freqtuenclee  between  the  cut¬ 
off  frequency  and  ^  vhen  considering  any  particular  mode.  With 
this  restriction,  we  may  speak  of  a  single-vedued  function 
for  each  t.  This  function  is  obtained  in  principle  by  evaluating 

As  before,  we  may  describe  this  function  by  two  parametric 

equations  giving  and  t  in  teims  of  ahese  equations  may  then 

be  used  for  numerical  congputatlons  or  as  a  basis  of  analytic 

studies.  To  demonstrate  that  numerical  computations  of  this 

function  are  feasible,  we  give  in  Fig.  l6  a  plot  of  ^^(t)  for 

2 

several  values  of  n,  for  D  >  1  and  A  >■  O.5625,  idilch  was  conq^uted 
In  this  msumer. 

The  qualitative  fom  of  Eq(t)  Is  clear  from  our  previous  dis¬ 
cussion  of  E_((u)  and  ui  _(t).  It  Is  zero  at  r  ■  1,  rises  to  a 
a  n^g 


maxlmvBi,  and  then  decreases.  Ihe  decrease  vlll  be  interrupted  by 

a  narrow  singularity  near  t  ..  However,  this  latter  feature 

n,r 

shall  be  ignored,  since  we  are  neglecting  the  Airy  phase. 

It  is  instructive  to  examine  Ej^(t)  in  the  limit  of  large  n. 

In  this  limit  we  may  derive  a  single  closed  expression  for  E^. 

We  find 

- Ho  eXP  {-  D(t-1)/[a2-(t-1)]}.  (6.ll^) 

^  Bg[A  -(t-I)]^/*^  ■* 

The  derivation  of  this  e^^resslon  is  given  in  Appendix  C. 

The  most  Interesting  feature  of  (6.l4)  Is  that  it  depends  on 
n  only  throu^  the  factor  B~^.  Thus  the  time  of  maximum  amplitude 
for  a  high  frequency  source  is  independent  of  mode  number.  We 
find  this  to  occur  \hen 

1/2 

T  -  1  =  A^Cl  +  ^]  +  3/4  -  {(|a^3/4)^  +  BA^}  .  (6.15) 

2 

If  A  is  small  (shallow  wave  guide)  or  D  is  large,  this  is 
approximately 

T  -  1  -  A^(1+(2/3)A^)‘^.  (6.16) 

In  general,  t  decreases  if  D  is  increased. 

We  may  estimate  the  duration  of  the  wave  by  setting  the  ex¬ 
ponent  in  (6.14)  equal  to  some  constant  and  solving  for  t.  If  we 
take  this  constant  to  be  10,  we  find  that  the  value  of  is 
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10(10+D)^^  -10 

- - 


(6.17) 


cuid 


(t-1) 


IQA 

lO+D 


(6.18) 


We  recall^  from  our  previous  discussion,  that  the  maximum  value  of 

E  was  of  the  order  of  (B  A)”^  for  D  “  1.  The  quantity  (6.17)  is 
n  c 

clearly  much  less  than  this  for  emy  appreciable  vedue  of  D.  Thus, 
(6.18)  affords  a  reasonable  estimate  of  the  duration  of  the  wave. 
Discussion 

The  general  method  we  have  used  for  studying  the  characteris¬ 
tics  of  the  normal  modes  Is  based  on  the  fact  that  we  may  formally 
describe  the  variation  of  any  characteristic  with  respect  to  fre¬ 
quency  or  T  by  giving  two  parametric  eq,\aatlons.  In  addition  to 
the  feasibility  of  the  method  for  numerical  calculations.  It  has  a 
distinct  advantage  as  a  starting  point  for  analytical  studies. 

One  may  qualitatively  describe  the  variation  of  the  characteristic 
by  examining  the  form  of  the  parametric  equations.  !Bius,  for 
example,  we  proved  that  the  group  velocity  had  only  one  mlnlmvan. 
When  we  asked  more  quantitative  questions,  the  study  of  the  para¬ 
metric  eqviations  suggested  means  of  deriving  approximate  euiswers 
to  these  questions.  Formerly,  some  of  these  could  only  be  answered 
after  a  lengthy  nisnerlual  calculation.  As  an  example,  we  derived 
an  expression  for  the  amplitude  of  the  Airy  phase  and  showed  that 
this  was  negligible  c(»ipared  to  the  maximuu  amplitude  In  the  ground 
wave  for  higher  order  modes. 


VII.  qUALUATIVE  DESCRIPTKW  OF  THE  WAVEFOBKS 

In  this  section,  we  shall  apply  the  results  of  the  preceding 
sections  to  the  Interpretation  of  the  vavefoms  which  would  be 
observed  at  relatively  large  dlstemces. 

Comparison  of  the  Contributions  from  the  Individual  Modes 

The  contribution  to  the  waveform  from  a  single  normal  mode  Is 
a  wave  of  slowly  varying  freq,uency.  As  we  shoved  In  the  previous 
section,  the  wave  Is  generally  composed  only  of  a  ground  wave.  In 
this  event,  we  may  associate  a  single  frequency  with  each  value 
of  T  (time  in  units  of  c^r).  'Ihis  was  a  direct  result  of  placing 
source  and  receiver  outside  the  guide.  !Rxe  amplitude  of  the  Airy 
ihase  may  be  appreciable  for  the  lowest  order  modes,  but  this 
amplitude  decreases  rapidly  If  one  Increases  D. 

Ihe  principal  effect  of  source  and  receiver  being  outside  the 
guide  was  to  restrict  the  band  width  carried  by  any  given  mode, 
nils  Is  clear  from  our  studies  of  and  E^Cui)  In  the  previous 

section.  Bie  exponential  factor  e'*®  ^  la  the  relative 

excitation  amplitude  decreased  rapidly  for  frequencies  above  the 
cutoff  frequency  for  almost  any  appreciable  value  of  D.  The  band 
width.  In  the  limit  of  large  mode  nunber,  was  very  small,  and 
proportional  to  n”^.  In  this  limit  there  wovild  be  very  little 
overlap  of  the  bands  carried  by  different  modes.  However,  seme 
overlap  would  occur  for  the  lowest  order  modes  If  D  were  not  large. 

The  contribution  from  a  given  mode  will  depend  on  the  fre¬ 
quency  spectrvBn,  g(w},  of  the  source.  However,  If  g(u))  does  not 
vary  significantly  over  the  band  of  frequencies  carried  by  the 
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mode,  it  Is  clear  that  one  may  tcdte  g((u)  as  being  constant  over 
this  fre(iuancy  band.  Ihls  will  be  true  for  higher  order  modes 
especially.  Qien  one  may  assume  g(tt>)  to  be  equal  to  Its  value  at 
the  cutoff  frequency  as  a  reasonable  approximation.  For  this  to 
be  valid  for  lowest  order  modes,  we  should  assvime  that  g(tt))  does 
not  vaiy  significantly  over  frequency  ranges  of  the  order  of  uu. 

This  would  be  true  for  any  realistic  time  dependence  of  the  source, 
provided  that  the  signal  duration  near  the  source  were  less  than 
l/tt).  For  example.  If  f(t)  Is  sin  art  between  t  »  0  and  t  -  2n/cr, 
and  otherwise  zero,  this  criterion  will  certainly  be  true  If  or 
Is  of  the  order  of  10  w. 

The  time  duration  of  the  contribution  from  a  given  mode,  at 
a  given  distance  r,  will  be  roughly  Independent  of  mode  number. 

This  Is  clear  from  equation  (6.14),  which  shows  that,  in  the  limit 
of  large  n,  the  relative  excitation  amplitude  factors  Into 
euid  a  quantity  that  depends  only  on  t.  This  will  not  be  true  for 
the  low  order  modes,  as  the  signal  duration  will  then  be  limited  by 
the  parameter  which  increases  with  mode  number. 

It  follows,  from  (6.14)  and  our  previous  remarks  concerning 
g(ui),  that  the  amplitudes  of  the  various  normal  mode  waves,  at 
any  given  time,  will  vary  with  n  as  B^^g(u)^),  for  sufficiently 
large  n. 

Dominant  Frequencies 

If  one  were  to  plot  the  absolute  values  of  the  wave  amplitudes 
for  the  various  modes,  corresponding  to  a  given  t,  versus  mode 
nunber,  he  would  find  that  the  resulting  curve  has  a  maximum  for 
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Bone  value  of  n.  Ibis  value  of  n  would  give  the  domlnwt  mode 
corresponding  to  that  particular  t.  IRie  frequency  would  be 

the  dominant  frequency  at  that  given  value  of  t.  This  frequency 
will  clearly  not  be  that  frequency  for  which  the  absolute  veuLue  of 
g((u)  Is  a  maximum. 

In  the  event  that  g(u))  has  Its  maximum  at  a  frequency  large 
In  comparison  to  Wf  and  varies  relatively  smoothly  with  frequency, 
our  previous  remarks  Indicate  that  this  domlmmt  frequency  will  be 
less  than  that  for  idilch  g(iu)  Is  maximum.  As  an  example,  we  con¬ 
sider  the  hypothetical  case 

|g((u)|  »  .  (7.1) 

The  maximum  value  occurs  when  w  *  Q.  We  assume  Q  »  u).  It  Is 
clear  that  emy  source  with  such  a  spectrum  will  primarily  excite 
modes  of  high  order.  Ihus,  one  may  use  (6.14)  euid  approximately 
take  the  amplitude  of  the  n-th  normal  mode  as  being  proportional 
to  B^^|g(u))|  or  ,  where  u)^  is  the  cutoff  frequency  of 

the  n-th  normal  mode.  Ihe  dominant  normal  mode  Is  then  that  whose 
cutoff  frequency  Is  approximately  Ihe  dominant  frequency 

will  be  of  this  order  of  magnitude  and  hence  slgnlflccuitly  less 
than  q,  the  dominant  frequency  In  the  source  spectrum. 

It  follows  that  the  signal  observed  at  large  distances  will 
give  the  appearance  of  being  ccmq>osed  of  frequencies  from  the  low 
frequency  spectrum  of  the  source.  The  higher  frequencies  will 
appear  to  have  been  "eliminated. " 
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!Ehe  above  renarka  afford  us  vith  a  certain  qLualltatlve  inslflkt 
Into  the  gross  features  of  the  vavefoxna.  Haturallyi  no  single 
mode  will  doalnate,  but,  rather,  a  group  of  nodes  vlll  be  more 
predominant  than  others* 

Bie  Bgplrlcal  Rxase  Velocity 

Ihe  definition  of  phase  velocity,  as  usually  given,  Is 
applicable  only  to  a  single  normal  mode.  However,  one  may  give  an 
K^irlcal  definition  of  phase  velocity  as  the  apparent  <veloclty 
with  \dilch  characteristic  features  in  the  waveform  move  In  the 
horizontal  direction.  This  definition  Is  often  applle<l  In  experi¬ 
ments  on  anomalous  sotmd  propagation  In  the  atmosphere  to  determine 
the  angle  of  Incidence  of  the  "ray"  to  the  ground. This 
definition  vlll  be  meaningless,  however,  for  long  distance  propa¬ 
gation,  unless  the  distortion  of  the  waveform  with  distance  Is 
small.  This  Is  clear  since  the  waveform  vlll  have  a  long  duration 
and  it  may  be  difficult  to  decide  which  peak  In  the  waveform  observed 
at  one  distance  corresponds  to  a  certain  peak  in  the  waveform  ob¬ 
served  at  another  distance.  Iftiless  there  Is  a  reasonable  amount 
of  correlation  between  waveforms  observed  at  different  distances, 
this  question  may  be  virtually  Impossible  to  answer. 

Since  each  mode  has  a  different  variation  of  phase  velocity 
with  T,  one  would  expect.  In  the  general  case,  when  many  modes 
contribute  to  the  waveform,  that  the  correlation  between  waveforms 
observed  at  different  distances  would  be  small,  euad  that  the 
enqplrlcal  phase  velocity  would  be  111  defined.  For  the  particular 
model  considered  in  this  paper,  with  source  and  receiver  outside 


-58- 


tbe  guide,  the  distortion  of  waveform  with  distance  Is  relatively 
small,  nie  reason  for  this  clrcianstance  Is  clear  when  one  considers 
how  the  phase  velocity  varies  with  frequency.  For  any  mode,  the 
phase  velocity  is  c^  at  the  cutoff  frequency  and  decreases  relatively 
slowly  (as  contrasted  to  the  ero\xp  velocity)  with  increasing  fre¬ 
quency.  Since,  with  source  and  receiver  outside  the  guides,  each 
mode  carries  only  frequencies  close  to  Its  cutoff  frequency,  the 
phase  velocity  will  be  very  nearly  constant  for  all  frequencies  of 
appreciable  amplitude. 

For  each  value  of  t,  there  is  an  inherent  uncertainty  in  the 

empirical  phase  velocity.  If  the  source  excites  the  modes  for  n 

between  n^  and  n^,  the  empirical  phase  velocity  could  range  from 

V  (t)  to  V  (t),  depending  on  Just  what  cJiaracteristic  feature 
“l 

in  the  waveform  is  being  used  as  a  basis  for  the  measurement.  The 

difference  v  -  v  may  be  considered  as  a  velocity  of  distortion. 

^1  ^2 

If  observations  are  made  at  r  and  im-L,  the  correlation  between  these 

wavefozms  at  corresponding  times  may  be  expected  to  be  small  if  L 

is  greater  or  of  the  order  of  magnitude  of  v  u)  (v^  ”  \  ^  >  where 

(u  is  the  dominant  frequency  at  the  corresponding  t,  and  v  is  the 

average  phase  velocity  of  all  the  excited  modes  at  the  same  t.  If 

L  is  this  large,  the  phase  difference  between  the  n^^-th  and  n^-th 

normal  mode  waves  will  have  changed  by  roughly  one  radian. 

We  may  therefore  make  a  rough  estimate  of  the  distance  the  wave 

must  travel  before  the  waveform  is  significantly  distorted.  Assuming 

the  soiuve  primarily  excites  high  order  modes,  we  use  (5*33)  to 

estimate  the  difference  v  -  v^  .  Bie  critical  factor  in  the 

n,  n_ 
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resultlng  expression  will  be  All  other  factors  nay  be  re¬ 

placed  by  averages  or  order  of  magnitude  e<iulvalents  without  sig¬ 
nificantly  changing  L.  Thus  v^  and  are  replaced  by  c^,  n^^  +  Ug 
Is  r^laced  by  4u)/(TT(ij)^  where  u>  Is  the  dominant  frequency  of  the 
source,  etc.  Furthermore,  assmlng  that  a  large  number  of  modes 
are  excited,  we  nay  take  ng  -  n^  as  being  (2/tt)(A(u/u)),  where  Am  Is 
the  band  width  of  the  source  spectrum.  Doing  this,  we  find 


(1^^A^  -  T 
(t  -  1)2a 


H 


(7.2) 


As  would  be  expected,  the  length  L  Is  very  large  for  t  close 
to  1  and  decreases  with  Increasing  t.  It  Is  large  for  higher  fre¬ 
quency  sources  or  sources  of  narrow  band  width. 

Variation  of  Waveform  with  Time  at  a  Given  Distance 

Since  the  frequency  for  the  ground  wave  Increases  with  t 

for  each  normal  mode.  It  Is  clear  that  the  earliest  portions  of  the 
waveform  must  be  characterized  by  lower  frequencies  than  are  present 
In  the  later  portions  of  the  waveform.  It  Is  Interesting  to  rote, 
however,  that,  for  sources  ^Ich  primarily  excite  higher  order  modes, 
the  frequency  spectrum  of  the  waveform  Is  nearly  uniform  over  the 
duration  of  the  signal.  tRils  Is  clear  from  (6.l4),  which  shows  the 
variation  of  relative  excitation  amplitudes  with  t  to  be  Independent 
of  mode  number. 

Furthermore,  for  hl^er  order  modes,  the  frequency  changes  very 
slowly  with  T.  nils  Is  clear  from  (^.30).  Ihe  waveform  as  a  whole 
may  be  expected  to  have  meocimum  amplitudes  at  the  value  of  t  given 
by  (6.l6),  since  the  relative  amplitudes  of  all  the  higher  order  modes 


will  be  maxlmvBi  at  this  value*  She  variation  of  peak  an^litudes 
with  T  vill  be  roughly  as 


throughout  the  duration  of  the  signal.  This  also  follows  from  (6.14). 

If  lover  order  modes  are  excited,  their  effect  on  the  total 
vavefom  will  be  most  noticeable  during  the  earlier  portion  of  the 
signal.  Ihls  is  clear  from  Pig.  l6,  which  shows  that  the  lover 
order  modes  have  their  peak  amplitudes  at  earlier  times  than  the 
limit  represented  by  (6.I6).  At  later  times,  these  modes  will  not 
be  present  at  «J.l,  since  their  duration  is  limited  by  which 

is  significantly  less  than  the  maximum  value  of  t  at  which  the 
signal  may  be  observed. 
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vm.  sutmcf  and  coHcmpiiio  aamacs 

Ve  have  applied  the  method  of  noznal  modes  to  a  slsqple  model 
of  layered  media,  for  the  particular  case  when  both  source  and  re¬ 
ceiver  were  outside  the  layer  with  lowest  sound  velocity.  Ihe  diffi¬ 
culties  usually  encountered  in  applying  the  method  In  guided  wave 
problems  were  alleviated  by  the  Introduction  of  a  variable  Oils 
enabled  us  to  write  the  dispersion  relations  in  a  parametric  form. 
Furthermore,  It  enabled  us  to  answer  some  q,uestlons  of  theoretical 
Interest,  such  as  that  of  the  location  of  the  complex  roots  In  the  k- 
plane. 

It  was  found  that  the  variation  of  Euiy  characteristic  of  the 
normal  modes  with  frequency  could  be  described  by  two  parametric 
equations  In  terms  of  u*  parametric  equations  were  found  to  be 
readily  amenable  to  analytic  investigation  and  we  were  able  to  study 
the  variation  of  these  characteristics  with  mode  nvnber.  Gb  the  basis 
of  this  study,  ve  were  able  to  make  some  qualitative  predictions  of 
the  nature  of  the  waveforms  which  would  be  observed  at  large  dis¬ 
tances  (Sec.  VIl)  when  a  large  number  of  modes  were  excited. 

Ihe  methods  ve  have  used  for  studying  the  characteristics  of 
the  normal  modes  will.  In  general,  be  Inapplicable  to  more  complex 
models  of  layered  media.  It  Is  hoped  however,  that  the  results  ob¬ 
tained  for  this  model  may  be  of  use  In  the  development  of  a  deeper 
Insight  Into  the  nature  of  guided  waves. 
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Holfspoct  Co 


Diagram  of  the  model  considered  In  this  paper. 


Fig.  3  -  Allowed  values  of  in  veurlous  regions  of  the 

k-plane.  (it  is  assumed  that  -n  <  <  )t>) 

■“  R  “■ 


Fig.  5  -  Plot  of  path  traced  by  roots  corresponding  to 


given  m  in  liL^c^/u  -  plane  as  u  is  varied. 

(A^»  0.25  or  c,/c  =  0.9.) 
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Frequency  in  units  of  bJ 

Pig.  7  -  Plot  of  (units  of  2/h)  vs.  u  (units  of  u)  for  the 

complex  nxxies.  (A^  »  0.2^.}  The  curves  are  shovm  only 
for  frequenles  between  the  two  cut-off  frequencies  of  the 
corresponding  mode. 
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Frequency  in  units  of  S 


Fig.  8  -  Sketch  of  group  velocity  vs.  B  =  w/u  showing 
characteristic  behavior. 


Mode  number  n 


Fig.  9  -  Plot  of  6b  ,  vs.  n  for  c,/c  =  0.8. 
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Modt  numbtr  n 


Flg«  11  -  Freqpieacy  relative  to  cutoff  frequency  at  group 
velocity  mlninum  vs.  mode  number  n.  Frequency  is 
in  units  of  u.  =  0.8.) 
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Fig.  15  -  Sketch  of  the  variation  of  d  k^^/du)  with  frequency 
f or  n  =  1  and  n  >  1. 
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Frequency  (o  (units  of  IS) 


Fig.  l4  -  Sketch  of  the  variation  of  Gjj(w)  with  B  =  w/<*> 
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Fig.  15  -  Sketch  of  the  variation  of  E°(u))  with  frequency. 
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APPENDIX  A 


Equation  (3*^)  Is  satisfied  if  and  only  if 
B  cos  g,  -  M,  =  ^  (m  +  l)Tr 


(A.l) 


for  some  integer  m.  Ibis  is  equivalent  to  two  simultaneous  equations 
if  IX  is  complex.  Ihese  two  equations  may  be  taken  in  the  form: 


sin  I 


tan  5 


B  sinh  J\* 


_ i _ !L_ 

[i(m.l)n  -  51 


(A.  2) 


(A.3) 


i 


I 

I 


xbere  we  have  set  |i  =  ■  5  *  IT)* 

The  first  equation  may  be  solved  eaqjlicltly  for  I  if  the  right 

hand  side  is  less  than  1.  (if  the  right  hand  side  is  greater  ttuui 

1,  no  solution  will  exist.)  One  may  therefore  use  (A. 2)  and  (3<8) 

2  2 

to  e]q>ress  k  in  terms  of  7].  Since  k  is  complex,  we  will  have  two 

2  2 

parametric  equations  giving  both  (k  )p  and  (k  in  terms  of  7^. 

One  may  then  let  T|  run  over  a  sequence  of  values  and  compute  (& 

Q 

and  (k  for  each  1].  Iben  one  may  compute  the  corresponding  values 
of  kj^  and  k^  from  the  relations 

f  1  2  1  gS  p2 

.  {.  i  .  I  [(k=)^  .  (k®)j  ]  }  ,  (A.i.) 

1 

\ 

where  the  plus  sign  is  used  to  find  k^  and  the  minus  si^  to  find  | 

k^.  In  this  manner  the  line  (3*H)  °iay  be  tabulated. 


1 


Ibe  eq,uatlon  (A>3)  loay  be  easily  solved  aunerlcedly  for 
This  is  done  most  simply  by  taking  the  arc  tangent  of  both  sides 
and  then  iterating,  starting  with  any  cmvailent  value  of  §  between 
0  and  tt/2.  abe  convergence  will  be  rapid  for  m  >  1.  (One  may 
show  that  there  will  be  no  solution  of  (A* 3)  for  which  0  <  5  <  n/2 
unless  m  >  1.  The  iteration  method  is  tlierefore  satisfactory  for 
all  cases  of  interest.)  One  may  let  1)  run  through  a  seq.uence  of 
values  and  compute  ^  for  each  1).  This  value  of  §  is  then  used  to 
compute  B  «  u/S  from  (A. 2).  Then  (k  )g  and  (k  may  be  computed 
by  use  of  (3.6).  Finally,  k^  wd  k^  may  be  computed  by  use  of  (A.^)< 
For  each  value  of  7],  one  tabulates  w,  k^,  cmd  k^. 

The  curves  of  k^  vs.  kj,  or  kj^  vs.  lu,  or  k^  vs.  tu  may  then  be 
plotted.  This  was  the  method  used  to  find  the  curves  shown  in 
Figures  5>  6,  and  7* 

The  two  equations,  (3*12)  and,  (3*15)#  were  obtained  after  an 
examination  of  the  equations,  (A. 2),  (A. 3)^  (3>6),  and  (A. 4),  in 
the  limit  of  large  T). 


! 

I 
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AFPENDIX  B 

Bie  quantity  k^”  is  given  by  the  e:q>re8Bion 

K'  ■  t3  Tbb  *  B  Tjjjl,  (B.1) 

vhere 


^  ^BB 

(B.2) 

^BBB  “  ^4  ^BB  * 

(B.3) 

2 

=  A  (sin  u  cos  \i)/T, 

(B.4) 

T  »  CA^(1  -  2  sin^n)  -  T  ]/T, 
piki  '  '  a  '  ' 

(B.5) 

T  =  -  T  (4T  +  3T  )/T, 

M-UU  M-  UU 

(B.6) 

T  =  (1  +  A^  Bin^ji)^/^, 

(B.7) 

^ 

(B.8) 

(B.9) 

‘*BBB  “  ■  ®UM.h  '  ^^*3  ^BB 

B  , 

(B.IO) 

B^  ■  (l  +  B  sin  (i)/coB  \ii, 

(B.ll) 

B  -  B  +  2B  tan  u. 

(B.12) 
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B  -  3B  +  (3B  -  B)  tan  n.  (B.I3) 

To  tabulate  k^"  vs.  uj,  one  lets  u  run  through  a  sequence  of  values* 
For  each  value  of  he  computes  B  from  (3*9)  and  then  proceeds  to 
calculate  k^"  by  computing  the  relevant  quantltltes  In  the  following 
sequence:  T,  B^,  B^^^,  B^^^,  ^ig,  Uggg,  Tgg,  Tggg, 

k^"  .  This  scheme  could  probably  be  simplified,  but  doing  this 
would  require  extensive  algebraic  manipulations.  The  method  we 
outline  above  leaves  the  algebra  to  an  electronic  computer. 


AFFEXIDIX  C 


If  «  1, 


oae  nay  approximately  take 


c 

V 


o 


n 


T 


1  + 


B  g, 

1  +  B^g 


or,  eq,uivaleiitly. 


M- 


(C.l) 


(C.2) 


If  one  examines  the  expression  (5*2)  for  the  phase  velocity  In  the 
limit  of  small  g,  he  may  obtain  equation  (^>33)  using  (C.2). 

To  derive  (6.14),  one  first  expresses  Ej^(u))  In  terms  of  li, 

-1/3 

assuming  u  «  B  '  Qae  finds,  in  lowest  order,  that 

V 

•  1  l/P  ** 

\  (l+B^n)-^'*^  B^g  e  .  (C.3) 

Equation  (6.l4)  follows  with  the  substitution  (C.2). 
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